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Using 347.5 fb^'^ of data recorded by the BABAR detector at the PEP-II electron-positron collider, 
244 X 10^ signal events for the — s- -k^ e'^ Ve decay channel are analyzed. This decay mode 
is dominated by the K (892)"^ contribution. We determine the K (892)° parameters: mj^.^ggjjo = 
(895.4±0.2±0.2) MeV/c^, r^^.^ggj^o = (46.5±0.3±0.2) MeV/c^ and the Blatt-Weisskopf parameter 
rBW ~ 2.1±0.5±0.5 ( GeV/c)^^ where the first uncertainty comes from statistics and the second from 
systematic uncertainties. We also measure the parameters defining the corresponding hadronic form 
factors at = {rv = = 1.463±0.017±0.031, r2 = = 0.801±0.020±0.020) and the value 

of the axial- vector pole mass parameterizing the variation of and A2: tua = (2.63±0.10±0.13) 
GeV/c^ The S -wave fraction is equal to (5.79±0.16±0.15)%. Other signal components correspond 
to fractions below 1%. Using the — > K^-k^-k^ channel as a normalization, we measure the D"*" 
semileptonic branching fraction: B{D+ K'n+e+Ue) = (4.00 ± 0.03 ± 0.04 ± 0.09) x 10"^ where 
the third uncertainty comes from external inputs. We then obtain the value of the hadronic form 
factor Ai at = 0: Ai{0) = 0.6200 ± 0.0056 ± 0.0065 ± 0.0071. Fixing the P-wave parameters we 
measure the phase of the S-wave for several values of the Kn mass. These results confirm those 
obtained with K-w production at small momentum transfer in fixed target experiments. 

PACS numbers: 13.20.Fc, 12.38.Gc, ll.30.Er, ll.15.Ha, 14.40.Df 



I. INTRODUCTION 

Detailed study of the £>+ — )• /r~7r+e+i^e decay channel 
is of interest for three main reasons: 

• it allows measurements of the different Ktt reso- 



*Now at Temple University, Philadelphia, Pennsylvania 19122, 
USA 

tAlso with Universita di Perugia, Dipartimento di Fisica, Perugia, 

Italy 

t Also with Universita di Roma La Sapicnza, 1-00185 Roma, Italy 
^Now at University of South Alabama, Mobile, Alabama 36688, 
USA 

^Also with Universita di Sassari, Sassari, Italy 



nant and non-resonant amplitudes that contribute 
to this decay. In this respect, we have measured 
the S'-wavc contribution and searched for radially 
excited P-wave and D-wave components. Accu- 
rate measurements of the various contributions can 
serve as useful guidelines to _B-meson semileptonic 
decays where there are still missing exclusive final 
states with mass higher than the D* mass. 

• High statistics in this decay allows accurate mea- 
surements of the properties of the K (892)° me- 
son, the main contribution to the decay. Both reso- 
nance parameters and hadronic transition form fac- 
tors can be precisely measured. The latter can be 
compared with hadronic model expectations and 
Lattice QCD computations. 
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• Variation of the Ktt S'-wavc phase versus the Kir 
mass can be determined, and compared with other 
experimental determinations. 

Meson-meson interactions are basic processes in QCD 
that deserve accurate measurements. Unfortunately, me- 
son targets do not exist in nature and studies of these in- 
teractions usually require extrapolations to the physical 
region. 

In the Ktt system, S-wave interactions proceeding 
through isospin equal to 1/2 states are of particular in- 
terest because, contrary to exotic 1 = 3/2 final states, 
they depend on the presence of scalar resonances. Stud- 
ies of the candidate scalar meson k = (800) can thus 
benefit from more accurate measurements of the 1 = 1/2 
S'-wavc phase below niKTv = 1 GeV/c^ [1]. The phase 
variation of this amplitude with the Kn mass also enters 
in integrals which allow the determination of the strange 
quark mass in the QCD sum rule approach [2, 3]. 

Information on the Kn S-w&ve phase in the isospin 
states 7=1/2 and 7 = 3/2 originates from various ex- 
perimental situations, such as kaon scattering, D — > Kttt: 
Dalitz plot analyses, and semileptonic decays of charm 
mesons and r leptons. In kaon scattering fixed tar- 
get experiments [4, 5], measurements from LASS (Large 
Aperture Solenoid Spectrometer) [5] start at mnTr = 
0.825 GeV/c^, a value which is 0.192 GeV/c^ above 
threshold. Results from Ref. [4] start at 0.7 GeV/c^ 
but are less accurate. More recently, several high statis- 
tics 3-body Dalitz plot analyses of charm meson hadronic 
decays have become available [6-9] . They provide values 
starting at threshold and can complement results from K 
scattering, but in the overlap region, they obtain some- 
what different results. It is tempting to attribute these 
differences to the presence of an additional hadron in 
the final state. The first indication in this direction was 
obtained from the measurement of the phase difference 
between S- and P- waves versus rriKTr in B'^ — ^ J/ipK~ Tr~^ 
[10] which agrees with LASS results apart from a rela- 
tive sign between the two amplitudes. In this channel, 
the J/ip meson in the final state is not expected to inter- 
act with the Ktt system. 

In T decays into Kwi/r there is no additional hadron 
in the final state and only the 7=1/2 amplitude con- 
tributes. A study of the different partial waves requires 
separation of the r polarization components using, for 
instance, information from the decay of the other r lep- 
ton. No result is available yet on the phase of the Ktt 
S-wave [12] from these analyses. In 7)+ K^ir^e^Vf, 
there is also no additional hadron in the final state. All 
needed information to separate the different hadronic an- 
gular momentum components can be obtained through 
correlations between the leptonic and hadronic systems. 
This requires measurement of the complete dependence 
of the differential decay rate on the five-dimensional 
phase space. Because of limited statistics previous ex- 
periments [13-15] have measured an 5-wave component 
but were unable to study its properties as a function of 
the TCtt mass. We present the first semileptonic charm 



decay analysis which measures the phase of the 7 = 1/2 
Kit S-wave as a function of rriKv from threshold up to 
1.5 GeV/c^. 



TABLE I: Possible resonances contributing to Cabibbo- 
favored T)"*" semileptonic decays [16]. 



resonance 




B{X Ktt) 


mass 


width 


X 




% 


MeV/c^ 


MeV/c^ 


^o(800) (?) 


0+ 


100(?) 


672 ± 40 


550 ± 34 


it'*(892) 


1" 


100 


895.94 ± 0.22 


48.7 ±0.8 


Ji:i(1270) 


1+ 





1272 ± 7 


90 ±20 


7fi(1400) 


1+ 





1403 ± 7 


174 ± 13 


7i'*(1410) 


1" 


6.6 ±1.3 


1414 ± 15 


232 ± 21 


7i:o*(1430) 


0+ 


93 ±10 


1425 ± 50 


270 ± 80 


7^2* (1430) 


2+ 


49.9 ± 1.2 


1432.4 ± 1.3 


109 ±5 


7i:*(1680) 


1" 


38.7 ±2.5 


1717 ±27 


322 ± 110 



Table I lists strange particle resonances that can ap- 
pear in Cabibbo-favored semileptonic decays. = 
1+ states do not decay into Kn and cannot be observed 
in the present analysis. The 7ir*(1410) is a 1~ radial exci- 
tation and has a small branching fraction into Kn. The 
7sr*(1680) has a mass close to the kinematic limit and 
its production is disfavored by the available phase space. 
Above the 7i'*(892) one is thus left with possible con- 
tributions from the 7i'*(1430), 7i:*(1410) and 7C*(1430) 
which decay into Kn through S-, P- and D-waves, re- 
spectively. At low Kn mass values one also expects an 
S-wave contribution which can be resonant (k) or not. A 
question mark is placed after the k = Kq{800) as this 
state is not well established. 

This paper is organized in the following way. In Sec- 
tion II general aspects of the Kn system in the elas- 
tic regime, which are relevant to present measurements, 
are explained. In particular the Watson theorem, which 
allows the relating of the values of the hadronic phase 
measured in various processes, is introduced. In Sec- 
tion III, previous measurements of the S-wave Kn sys- 
tem are explained and compared. The differential decay 
distribution used to analyze the data is detailed in Sec- 
tion IV. In Section V a short description of the detector 
components which are important in this measurement is 
given. The selection of signal events, the background 
rejection, the tuning of the simulation and the fitting 
procedure are then considered in Section VI. Results 
of a fit which includes the S-wave and K (892)° sig- 
nal components are given in Section VII. Since the fit 
model with only S- and P-wave components does not 
seem to be adequate at large Kn mass, fit results for sig- 
nal models which comprise S+K (892)"±7C (1410)" and 
S + K* {892)° + K* {UlO f + D components are given in 
Section VIII. In the same section, fixing the parameters 
of the K (892)" component, measurements of the phase 
difference between S and P waves are obtained, for sev- 
eral values of the Kn mass. In Section IX, measurements 
of the studied semileptonic decay channel branching frac- 
tion, relative to the — >■ K~n^n^ channel, and of its 
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different components are obtained. This allows one to 
extract an absolute normalization for the hadronic form 
factors. Finally in Section X results obtained in this anal- 
ysis are summarized. 



TABLE II: Predicted values for scattering length and effec- 
tive range parameters. 

Parameter [18] [19] 

a,y^ (GcV^M 1.52 1.60 ±0.16 
b}/- ( GcA' 17.0 ■■il.2 ± 1.5 



II. THE Kn SYSTEM IN THE ELASTIC ^iTT- r 

REGIME REGION ( ) 5-59 7.0 ± 0.4 



The Kn scattering amplitude {Tk-k) has two isospin 
components denoted T^/^ and T^/^. Depending on the 
channel studied, measurements are sensitive to differ- 
ent linear combinations of these components. In — > 
K-n+e+Ue, T- K^^n-Pr and B° -S- J/V'iV'-7r+ de- 
cays, only the / = 1/2 component contributes. The 
7 = 3/2 component was measured in K~^p — >■ ii'+TT+n 
reactions [4] whereas K~p — )• K~'7T~^n depends on the 
two isospin amplitudes: Tk-„+ = 1{2T^/'^ + T^/"^). In 
Dalitz plot analyses of 3-body charm meson decays, the 
relative importance of the two components has to be de- 
termined from data. 

A given Kw scattering isospin amplitude can be ex- 
panded into partial waves: 

oo 

T\s,t,u) = 167r^(2£ l)P£(cos 6')4(s) (1) 

where the normalization is such that the differential K-k 
scattering cross-section is equal to: 

do^ ^ A \Ti{s,t,uf 

dl? s (167r)2 ' ^ 

where s, t and u are the Mandelstam variables, 6 is the 
scattering angle and P^(cos^) is the Legendre polynomial 

of order t. 

Close to threshold, the amplitudes t[{s) can be ex- 
pressed as Taylor series: 

Re4(5) = ^sTs {P*f (4 + H {p*f + O {p*f) , (3) 

where aj and bj are, respectively, the scattering length 
and the effective range parameters, p* is the K or i: mo- 
mentum in the Ktt center-of-mass (CM). This expansion 
is valid close to threshold for p* < niTr- Values of aj 
and 6^ are obtained from Chiral Perturbation Theory 
[17, 18]. In Table II these predictions are compared with 
a determination [19] of these quantities obtained from 
an analysis of experimental data on Kw scattering and 
TTTT KK. Constraints from analyticity and unitarity 
of the amplitude are used to obtain its behavior close 
to threshold. The similarity between predicted and fit- 
ted values of aj^^ and is a non-trivial test of Chiral 
Perturbation Theory [18]. 

The complex amplitude t[ (s) can be also expressed in 
terms of its magnitude and phase. If the process remains 
elastic, this gives: 

= - = ^-n^,^(.)e^^^^(«). (4) 



Using the expansion given in Eq. (3), close to the thresh- 
old the phase (5j(s) is expected to satisfy the following 
expression: 

5l{s)^{p*f+'{a + P{p*f). (5) 

Using Eq. (3), (4) and (5) one can relate a and /3 to a\ 
and h\: 

a = 4 and/3 = hi + '^{a[f5io. (6) 

In Eq. 6, the symbol 5io is the Kronecker 5 function: 

Joo = 1, 5io = for Z ^ 0. 

The Watson theorem [20] implies that, in this elastic 
regime, phases measured in Kn elastic scattering and in 
a decay channel in which the Kn system has no strong 
interaction with other hadrons are equal modulo n ra- 
dians [21] for the same values of isospin and angular 
momentum. In this analysis, this ambiguity is solved 
by determining the sign of the 5-wave amplitude from 
data. This theorem does not provide any constraint on 
the corresponding amplitude moduli. In particular, it is 
not legitimate (though nonetheless frequently done) to 
assume that the 5'-wave amplitude in a decay is propor- 
tional to the elastic amplitude t\{s). The Kn scattering 
S-wave, / = 1/2, remains elastic up to the Krj thresh- 
old, but since the coupling to this channel is weak [22], 
it is considered in practice to be elastic up to the Kr}' 
threshold. 

Even if the Kn system is studied without any accom- 
panying hadron, the S- or P-waves amplitudes cannot be 
measured in an absolute way. Phase measurements are 
obtained through interference between different waves. 
As a result, values quoted by an experiment for the phase 
of the 5-wave depend on the parameters used to deter- 
mine the P-wave. For the P-wave, the validity domain of 
the Watson theorem is a-priori more restricted because 
the coupling to Krj is no longer suppressed. However the 
p*^ dependence of the decay width implies that this con- 
tribution is an order of magnitude smaller than Kn for 
mK^< 1.2 GeV/c^. 

For pseudoscalar-meson elastic scattering at threshold 
all phases are expected to be equal to zero (see Eq. (5)). 
This is another important difference as compared with 
Dalitz plot analyses where arbitrary phases exist between 
the different contributing waves due to interaction with 
the spectator hadron. It is thus important to verify if 



8 



apart from a global constant S'-wavc phases measured 
versus rriKv: in 3-body D — )■ Ktttt Dalitz plot analyses, 
depend on the presence of the third hadron. Compar- 
ison between present measurements and those obtained 
in three-body Dalitz plot analyses are given in Section 
VIII B. 



III. PREVIOUS MEASUREMENTS 

In the following sections, we describe previous mea- 
surements of the phase and magnitude of the Kn S'-wavc 
amplitude obtained in K'^p scattering at small transfer, 
in T semileptonic decays, D meson three-body decays, 
and in charm semileptonic decays. 



A. Kn production at small momentum transfer 



A Kit partial wave analysis of high statistics data for 
the reactions K^p — >■ K^TT~^n and K^p — >■ K^n~ A'^~^ 
at 13 GeV, on events selected at small momentum trans- 
fer [4] , provided information on Kw scattering for rnxw in 
the range [0.7, 1.9] GeV/c^. The / = 3/2 Kn scattering 
was studied directly from the analyses of K^p — >• JsT+tt+ti 
and K~p — >■ i>sr~7r~A++ reactions. The phase of the 
elastic amplitude was measured and was used to 

extract the phase of the 1—1/2 amplitude from measure- 

1 /2 

ments of K^ir^ scattering. Values obtained for Sg are 
displayed in Fig. 1 for rriKTr < 1.3 GeV/c^, a mass range 
in which the interaction is expected to remain elastic. 
Above 1.46 GeV/c^ there were several solutions for the 
amplitude. 

A few years later, the LASS experiment analyzed data 
from 11 GeV/c kaon scattering on hydrogen: K~p 
K~'K^n [5]. They performed a partial wave analysis of 
1.5 X 10^ events which satisfied cuts to ensure Kn pro- 
duction dominated by pion exchange and no excitation 
of the target into baryon resonances. 

The Kn, I = 1/2, 5- wave was parameterized as the 
sum of a background term {BG) and the ifg (1430), which 
were combined such that the resulting amplitude satisfied 
unitarity: 



= sin 5 



1/2 JS)/^, 



BG ' 

+ e2'''BGsin5^.(i430) 



(7) 



g«'5jr*(i430) 



= Sin (^VJ + 5^5(,430))e<'-+'-S(-o)) ^ 

-1/2 

where and (5k* (1430) depended on the Kn mass. 

The mass dependence of was described by means 
of an effective range parameterization: 



cot {S]il) = 



1/2 

°'S,BG P 



+ 



l1/2 * 
bs,BG P 



(8) 



where a^^Q is the scattering length and b]/'^Q is the 
effective range. Note that these two parameters are dif- 
ferent from a| and bj introduced in Eq. (3) as the latter 
referred to the total amplitude and also because Eq. (8) 
corresponds to an expansion near threshold which dif- 
fers from Eq. (5). The mass dependence of Sk'{i43o) was 
obtained assuming that the Kq{1A30) decay amplitude 
obeys a Breit-Wigner distribution: 



cot ((5^.(1430)) 



m 



if* (1430) 



m 



if 71 



'^i£'*(1430)rif*(1430) (rriKn) ' 



(9) 



where 



*(1430) 



is the pole mass of the resonance and 



if* (1430) 



{juk-k) its mass-dependent total width. 



The total 7=1/2 S-wacve phase was then: 



.1/2 

^LASS 



,1/2 



BG +^K^(1430) 



(10) 



The LASS measurements were based on fits to mo- 
ments of angular distributions which depended on the 
interference between S-, P-, .waves. To obtain the 
I = 1/2 K~n~^ 5-wave amplitude, the measured I = 3/2 
component [4] was subtracted from the LASS measure- 
ment of TV-7r+ and the resulting values were fitted using 
Eq. (10). The corresponding results [23] are given in Ta- 
ble III and displayed in Fig. 1. 



TABLE III: Fit results to LASS data [23] for two mass in- 
tervals. 



Parameter 


rriKn € [0.825, 1.52] tuk 

GeV/c^ 


^ e [0.825, 1.60] 

GcV/c^ 


mK'iu-MY) (MeV/c^) 
rK*{i43o) (MeV/c^) 


1435 ± 5 


1415 ± 3 


279 ±6 


300 ±6 


4(sG (GcV-1) 


1.95 ±0.09 


2.07 ±0.10 


blilc (GeV-i) 


1.76 ±0.36 


3.32 ± 0.34 


B. 


T~ Km>r decays 





The BaBar and Belle collaborations [11, 12] measured 

the K'^n mass distribution in t~ — >■ K^n~Vr. Results 
from Belle were analyzed in Ref. [24] using, in addition 
to the ii'*(892): 

• a contribution from the iV'*(1410) to the vector 
form factor; 

• a scalar contribution, with a mass dependence 
compatible with LASS measurements but whose 
branching fraction was not provided. 

Another interpretation of these data was given in 
Ref. [25]. Using the value of the rate determined from 
Belle data, for the K*{lAlQi), its relative contribution to 
the D"*" K'n^e'^Ve channel was evaluated to be of 
the order of 0.5%. 
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FIG. 1: (color online) Comparison between the I — 1/2 S- 
wave phase measured in Kn production at small transfer for 
several values of the Ktt mass. Results from Ref. [4] are 
limited to rriKn < 1-3 GeV/c'^ to remain in the elastic regime 
where there is a single solution for the amplitude. The curve 
corresponds to the fit given in the second column of Table IIL 



C. Hadronic D meson decays 

Ktt interactions were studied in several Dalitz plot 
analyses of three-body D decays and we consider only 
D+ K-TT+TT+ as measured by the E791 [6], FOCUS 
[7, 8], and CLEO-c [9] collaborations. This final state is 
known to have a large iS-wave component because there 
is no resonant contribution to the ■k'^tt'^ system. In prac- 
tice each collaboration has developed various approaches 
and results are difiicult to compare. 

The S'-wave phase measured by these collaborations is 
compared in Fig. 2-a with the phase of the (/ = 1/2) 
amplitude determined from LASS data. Measurements 
from D'^ decays are shifted so that the phase is equal 
to zero for mK-K — 0.67 GeV/c^. The magnitude of the 
amplitude obtained in Dalitz plot analyses is compared 
in Fig. 2-b with the "naive" estimate given in Eq. (4), 
which is derived from the elastic (/ = 1/2) amplitude 
fitted to LASS data. 

By comparing results obtained by the three experi- 
ments analyzing _D+ — !■ i^^Tr+Tr+j several remarks are 
formulated. 



• A 7r+7r+ component is included only in the CLEO- 
c measurement and it corresponds to (15 ± 3)% of 
the decay rate. 

• The relative importance of / = 1/2 and / = 3/2 

components can be different in Ki: scattering and 



in a three-body decay. This is because, even if Wat- 
son's theorem is expected to be valid, it applies 
separately for the 7=1/2 and / — 3/2 compo- 
nents and concerns only the corresponding phases 
of these amplitudes. In E791 and CLEO-c they 
measured the total Ktt S-wave amplitude and com- 
pared their results with the / = 1/2 component 
from LASS. FOCUS [7], using the phase of the 
1 — 3/2 amplitude measured in scattering exper- 
iments, had fitted separately the two components 
and found large effects from the 1 — 3/2 part. In 
Fig. 2-a the phase of the total ^-wave amplitude 
which contains contributions from the two isospin 
components, as measured by FOCUS [8], is plotted. 

• Measured phases in Dalitz plot analyses have a 
global shift as compared to the scattering case (in 
which phases are expected to be zero at threshold) . 
Having corrected for this effect (with some arbi- 
trariness), the variation measured for the phase in 
three-body decays and in Kn scattering is roughly 
similar, but quantitative comparison is difficult. 
Differences between the two approaches as a func- 
tion of rriKn are much larger than the quoted un- 
certainties. They may arise from the comparison 
itself, which considers the total Ktt S-wave in one 
case and only the 7 = 1/2 component for scatter- 
ing. They could be due also to the interaction of 
the bachelor pion which invalidates the application 
of the Watson theorem. 

It is thus difficult to draw quantitative conclusions 



from results obtained with 



K- 



decays. 



Qualitatively, one can say that the phase of the S'-wave 
component depends roughly similarly on rn^^ as the 
phase measured by LASS. Below the K^{1430), the S'- 
wave amplitude magnitude has a smooth variation versus 
rriKTT- At the K^{\A3Q) average mass value and above, 
this magnitude has a sharp decrease with the mass. 



D. Di4, decays 

The dominant hadronic contribution in the 15^4 de- 
cay channel comes from the (J^ =1~) -ft'*(892) resonant 
state. E687 [13] gave the first suggestion for an additional 
component. FOCUS [14], a few years later, measured the 
5-wave contribution from the asymmetry in the angular 
distribution of the K in the Ktt rest frame. They con- 
cluded that the phase difference between S- and P- waves 
was compatible with a constant equal to 7r/4, over the 
-fi'*(892) mass region. 

In the second publication [26] they found that the 
asymmetry could be explained if they used the variation 
of the S'-wave component versus the Ktt mass measured 
by the LASS collaboration [5]. They did not fit to their 
data the two parameters that governed this phase varia- 
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rCeV/c') JGeV/c') 



FIG. 2: (color online) a) Comparison between the S-wave phase measured in various experiments analyzing the — s- K^tx^-k^ 
channel (E791 [6], FOCUS [7, 8] and CLEO [9]) and a fit to LASS data (continuous line). The dashed hue corresponds to the 
extrapolation of the fitted curve. Phase measurements from decays are shifted to be equal to zero at mK-n ~ 0.67 GeV/c^. 
b) The S'-wave amplitude magnitude measured in various experiments is compared with the elastic expression. Normalization 
is arbitrary between the various distributions. 



tion but took LASS results: 

1 



cot [6 



BG) — 



'^S.BG P 



(11) 



as^BG P* 2 
as.BG = (4.03 ±1.72 ±0.06) GeY~\ 
bs,BG = (1.29 ±0.63 ±0.67) GeV"^ 

These values corresponded to the total S'-wave ampli- 
tude measured by LASS which was the sum of / = 1/2 
and 7 = 3/2 contributions whereas only the former com- 
ponent was present in charm semileptonic decays. For the 
S'-wave amplitude they assumed that it was proportional 
to the elastic amplitude (see Eq. (4)). For the P-wave, 
they used a relativistic Breit-Wigner with mass depen- 
dent width [27]. They fitted the values of the pole mass, 
the width and the Blatt-Weisskopf damping parameter 
for the K*{S^2). These values from FOCUS are given in 
Table IV and compared with present world averages [16]. 
, dominated by the P-wave measurements from LASS. 

They also compared the measured angular asymme- 
try of the K in the Kn rest frame versus the K-k mass 
with expectations from a k resonance and conclude that 
the presence of a k could be neglected. They used 
a Breit-Wigner distribution for the k amplitude using 
values measured by the E791 collaboration [28] for the 
mass and width of this resonance {m^ = 797 ± 19 ± 
43 MeV/c2, r« = 410±43±87 MeV/c^). This approach 
to search for a k does not seem to be appropriate. Adding 



TABLE IV: Parameters of the /C*(892)° measured by FO- 
CUS are compared with world average or previous values. 

Parameter FOCUS results [26] previous results 
m^.o (MeV/c^) 895.41 ± 0.32l° g 895.94 ± 0.22 [16] 



r^f.o (MeV/c^ 
rsw (GcV/c)- 



47.79 ± 0.86 
3.96 ± 0.54: 



-1.32 
-1.06 
-1.31 
-0.90 



48.7 ±0.8 [16] 
3.40 ± 0.67 [5] 



a K in this way violates the Watson theorem as the phase 
of the fitted Kit amplitude would differ greatly from the 
one measured by LASS. In addition, the interpretation 
of LASS measurements in Ref. [19] concluded there was 
evidence for a k. In addition to the i^*(892) they mea- 
sured the rate for the non-resonant S-wave contribution 
and placed limits on other components (Table V). 

Analyzing D'^ — > K~TT'^e^Vf. events from a sample 
corresponding to 281 pb^^ integrated luminosity, the 
CLEO-c collaboration had confirmed the FOCUS result 
for the S'-wave contribution. They did not provide an 
independent measurement of the S'-wave phase [15]. 
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TABLE V: Measured fraction of the non-resonant S-wave 
component and limits on contributions from A'o(1430) and 
if* (1680) in the decay D+ K'-k+^j+v^,, obtained by FO- 
CUS [261. 



Channel 


FOCUS [26] (%) 


r(D + ->K-ir+M+i^j,) 


5.30 ± 0.74lg:^^ 


r{I5 + ^X-,r + M+i.,0;f*(1430) 


< 0.64% at 90% C.L. 


r(D + ->K-ir+M+i^^) 




< 4.0% at 90% C.L. 


T(D + ^K- ■K+ tJ.+ Uf^) 



a 



w 



V V| 




\ "^A 

V / 

Ve 1 


T% ^ ^ " 

\ \/ 



IV. 



K-Ti+e+Ve. DECAY RATE FORMALISM 



The invariant matrix element for the — > 
K~TT^e^Ue semileptonic decay is the product of a 
hadronic and a leptonic current. 



■n(jp^+)K{pK-)\s-i^ll - lb)c\D{pD+)) 

(12) 



V2 

Xw(Pi'e)7M(l - lh)v{Ve+) 



In this expression, px- , P7r+ 1 Pe+ ^-nd p^^ are the 
^ 7r+, and four- momenta, respectively. 
The leptonic current corresponds to the virtual 
which decays into e^v^. The matrix element of the 
hadronic current can be written in terms of four form 
factors, but neglecting the electron mass, only three are 
contributing to the decay rate: h and w±. Using the con- 
ventions of Ref. [29] , the vector and axial- vector compo- 
nents are, respectively: 



{'^{V^+)K{pK-)\sitiC\D{pD+)) 

= hef,ai3'rPD+ {Pk- + P^+f (Pk- - Pn+V ■ 

{■^{Pn+)K(j>K-)\sipi{-lb)c\D(j>D+)) 

= iw+ [pk- +P7t+)^ + iw^ [pk- - P^+ )f, ■ 



(13) 
(14) 



As there are 4 particles in the final state, the differ- 
ential decay rate has five degrees of freedom that can be 
expressed in the following variables [30, 31]: 

• TO^, the mass squared of the Kit system; 

• q^, the mass squared of the e'^Ve system; 

• cos (0k), where 6k is the angle between the K 
three-momentum in the Kw rest frame and the line 
of flight of the Kit in the D rest frame; 

• cos (Oe), where 9^ is the angle between the charged 
lepton three-momentum in the eVe rest frame and 
the line of flight of the ei^e in the D rest frame; 

• X: the angle between the normals to the planes de- 
fined in the D rest frame by the Kn pair and the 
eVe pair, x is defined between — tt and +Tr. 



FIG. 3: (color online) Definition of angular variables. 



The angular variables are shown in Fig. 3, where K^- 
is the K~ three-momentum in the Kit CM and is 
the three-momentum of the positron in the virtual W 
CM. Let i) be the unit vector along the Ktt direction in 
the D rest frame, c the unit vector along the projection 
of K K- perpendicular to ii, and d the vmit vector along 
the projection of perpendicular to ii. We have: 



m"^ = (P7r+ +PK-f', 9^ = (Pe+ +p^J^, (15) 



cos [Ok) = 



i) K 



K- 



v-K, 



cos ( 



Kk-\ ' ^ ^' \K,^\ ' 
cos (x) = c • d; sin (x) = (c x -D) • d 

The definition of x is the same as proposed initially 
in Ref. [30]. When analyzing D~ decays, the sign of x 
has to be changed. This is because, if CP invariance 
is assumed with the adopted definitions, x changes sign 
through CP transformation of the final state [14]. 

For the differential decay partial width, we use the 
formalism given in Ref. [29], which generalizes to five 
variables the decay rate given in Ref. [32] in terms of 
q^, cos 9k, cosOe and x variables. In addition, it pro- 
vides a partial wave decomposition for the hadronic sys- 
tem. Any dependence on the lepton mass is neglected as 
only electrons or positrons are used in this analysis: 



■X/3I{m^,q^ 



7K,0e 



,x) 



dm^dg^d cos {Ok)^ cos {9e)dx- 



(16) 



In this expression, X — pK-n tit. d where pktt is the mo- 
mentum of the Kt: system in the D rest frame, and 
P — 2p* /m. p* is the breakup momentum of the K-k 
system in its rest frame. The form factors h and w±, 
introduced in Eq. (13-14), are functions of m^, q^ and 
cos 9k ■ In place of these form factors and to simplify the 
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notations, the quantities J^i,2,3 are defined [29]: 



factors i?o,d 



T\ = X w+ + [;3 (px - + ) (Pe+ + Vv^ ) cos 



+ 



m 



if 



Ti = Pqmw-, (17) 
= pXqmh. 

The dependence of I on and x is given by: 

I = Ii +X2 cos 26*6 +X3sin2 6*6 cos 2x (18) 
+I4 sin 26e cos X + ^5 sin cos x 
+l6 cos + sin sin X 
+l8 sin 2^e sin x + sin^ sin 2x 

where Xi,...,9 depend on m^, q'^ and These quantities 
can be expressed in terms of the three form factors, J^i,2,3- 



r + ^in^^^dJ-^l^ 



•^31') I (19) 

1 



14 = (-^i-^^) sin^if 

15 = Re(J-rJ"3)sin0^f 

If = Im (J^iJ^2 ) sin^x 
Is = ^Ini(J"iJ"3)sin6'K 

Jg = -^Im(J-2J-3*)sin2^K 

Form factors -7^i,2,3 can be expanded into partial waves 
to show their cxphcit dependence on Ok- If only S-, P- 
and D- waves are kept, this gives: 



^1 = -^10 + ^11 cos 0K + T\ 



Scos^ 61k - 1 



1 

71 



^21 



Til cos 



(20) 



^3 = "^-^31 + Y ^^32 COS Qk ■ 

Form factors Tij depend on rr? and q^ . Txq characterizes 
the S'-wavc contribution whereas Ti\ and correspond 
to the P- and D-wave, respectively. 



A. P-wave form factors 

By comparing expressions given in Rcf. [29] and [32] it 
is possible to relate Ti\^ « = 1, 2, 3 with the helicity form 



Tx\ = 2V2aqHo 
T21 = 2aq{H++H_) 
7-31 = 2aq{H+-H_) 



(21) 



where a is a constant factor, its value is given in Eq. (26); 
it depends on the definition adopted for the mass distri- 
bution. The helicity amplitudes can in turn be related 
to the two axial- vector form factors Ai^iiQ^), and to the 
vector form factor V{q^): 



H±{q') 



1 



2m q 

mo + rn 



q^) {mo + m) Ax{q^) 

(22) 



[niD + m) Ai{q^) =F ; V{q^). 

mo +m 



As we are considering resonances which have an ex- 
tended mass distribution, form factors can also have a 
mass dependence. We have assumed that the and m 
dependence can be factorized: 

(V,Ai,^2)(9',m) = {V,AuA2){q'') x A{m) (23) 

where in case of a resonance A{m) is assumed to behave 
according to a Breit-Wigner distribution. 

This factorized expression can be justified by the fact 
that the q^ dependence of the form factors is expected to 
be determined by the singularities which are nearest to 
the physical region: q^ G [0, q^axl- These singularities 
are poles or cuts situated at (or above) hadron masses 
Mh — 2.1-2.5 GeV/c^, depending on the form factor. 
Because the variation range is limited to ~ 1 GeV^, 
the proposed approach is equivalent to an expansion in 
q^/M^j < 0.2. 

For the q"^ dependence we use a single pole parameter- 
ization and try to determine the effective pole mass. 



V{q') 

Ml') 

A2{q^) 



1-4- 

Ai(0) 



^2(0) 

1-4 



(24) 



where niv and ruA are expected to be close to mu* ~ 
2.1 GeV/c^ and mo^i — 2.5 GeV/c^ respectively. Other 
parameterizations involving a double pole in V have been 
proposed [33] , but as the present analysis is not sensitive 
to my, the single pole ansatz is adequate. 

Ratios of these form factors, evaluated at q"^ = 0, 
ry = A^(^Q-j and r2 = ^^|q| , are measured by study- 
ing the variation of the differential decay rate versus the 
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kinematic variables. The value of ^i(O) is determined by 

measuring the D"*" — K e+^e branching fraction. For 
the mass dependence, in case of the /i'*(892), we use a 
Breit-Wigner distribution: 



if* (892) 



?"lf(892)r°f.(892)Fi(m) 



'™l'*(892) -m'^ - im K' (892)^ K '{892) 



In this expression: 

• if* (892) is the K*{892) pole mass; 



(25) 



• ^^-(892) is tiie total width of the K*{892) for m = 

"lif*(892); 

• rx*(8g2)(m) is the mass-dependent K*{892) width: 

rif*(892)M = r°,.(3,2)S^^^J^i^(H; 

• Fi{m) = where B is the Blatt-Weisskopf 

damping factor: B = l/y^l + r^^yp*"^, rsw is the 
barrier factor, p* and Pq are evaluated at the mass 
m and mx*(892) respectively and depend also on 
the masses of the K*{892) decay products. 

With the definition of the mass distribution given in 
Eq. (25), the parameter a entering in Eq. (21) is equal 
to: 



a = 



3ttBk' 



(26) 



if* (892) 

where Bk' = B {K*{892) K-w+) = 2/3. 

B. S-wave form factor 

In a similar way as for the P-wave, we need to have 
the correspondence between the S-wave amplitude J^io 
(Eq. (21)) and the corresponding invariant form factor. 
In an 5-wave, only the helicity Hq form factor can con- 
tribute and we take: 



^10 = PKninD - 



1 - 



-As{m). 



(27) 



The term J^io is proportional to pkx to ensure that the 
corresponding decay rate varies as as expected from 
the L = 1 angular momentum between the virtual W and 
the 5'-wave Kw hadronic state. Because the variation 
of the form factor is expected to be determined by the 
contribution of = 1+ cs states, we use the same de- 
pendence as for Ai and A2. The term Asim) corresponds 
to the mass dependent 5- wave amplitude. Considering 
that previous charm Dalitz plot analyses have measured 
an 5-wave amplitude magnitude which is essentially con- 
stant up to the ^0(1430) mass and then drops sharply 



above this value, we have used the following ansatz: 

As = rsP{m)e'^''^"'\ and 
As = r-sP(mK5 (1430)) 

^ / ("^if* (1430) r^o* (1430))^ 

V (™^o*(1430) - '^^)' + (™ifo*(1430)rifJ(1430))' 



respectively for m below and above the i^o(1430) pole 
mass value. In these expressions, Ss{m) is the S'-wave 

phase, P{m) = 1 -|- r^g^ x a; -I- •* x + ... and x = 

— ^7 — )2 — 1. The coefficients r^Q have no dimension 



and their values are fitted, but in practice, the fit to data 
is sensitive only to the linear term. We have introduced 
the constant rs which measures the magnitude of the S'- 
wave amplitude. From the observed asymmetry of the 
cos 9k distribution in our data, rs < 0. This relative 
sign between S and P waves agrees with the FOCUS 
measurement [14]. 



C. D-wave form factors 

Expressions for the form factors J^,,2 for the I?-wave 
are [34]: 



^12 
^22 

^32 



moPKTT 



[{ml, -n? - q^) {mo + m) Ti{q^) 



T2{q^) 



^nriD mqpK^ {mo + m) Ti (g^). 



(29) 



2 2m\,mqp\^ 

3 {mD+fn) 



These expressions are multiplied by a relativistic Breit- 
Wigner amplitude which corresponds to the i^|(1430): 



^if„* = 



Td m^j(i430)r°f.(i430)^2(m) 



m 



if 1(1430) 



- im^.( 1430) r|fj( 1430) (m) 



.(30) 



ro measures the magnitude of the D-wave amplitude and 
similar conventions as in Eq. (25) are used for the other 
variables apart from the Blatt-Weisskopf term which is 
equal to: 



B2 = l/\l{rl^P*^-^f + 9rl^p*\ 
and enters into 



F2{m) 



B2{P*) 

PoJ MPoY 



(31) 



(32) 



The form factors Ti{q^) [i = 1, 2, V) are parame- 
terized assuming the single pole model with correspond- 
ing axial or vector poles. Values for these pole masses 
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arc assumed to be the same as those considered be- 
fore for the S- or P-wave hadronic form factors. Ratios 
of D-wAvc hadronic form factors evaluated at = 0, 
r22 = T2(0)/Ti(0) and r2v = Tv{0)/Ti{0) are supposed 
to be equal to one [35]. 

V. THE BABAR DETECTOR AND DATASET 

A detailed description of the BaBAR. detector and of the 
algorithms used for charged and neutral particle recon- 
struction and identification is providcid cilscwhere [36, 37]. 
Charged particles are reconstructed by matching hits in 
the five-layer double-sided silicon vertex tracker (SVT) 
with track elements in the 40 layer drift chamber (DCH), 
which is filled with a gas mixture of helium and isobu- 
tane. Slow particles which due to bending in the 1.5 T 
magnetic field do not have enough hits in the DCH, are 
reconstructed in the SVT only. Charged hadron identifi- 
cation is performed combining the measurements of the 
energy deposition in the SVT and in the DCH with the 
information from the Cherenkov detector (DIRC). Pho- 
tons are detected and measured in the CsI(Tl) electro- 
magnetic calorimeter (EMC). Electrons are identified by 
the ratio of the track momentum to the associated en- 
ergy deposited in the EMC, the transverse profile of the 
shower, the energy loss in the DCH, and the Cherenkov 
angle in the DIRC. Muons are identified in the instru- 
mented flux return, composed of resistive plate chambers 
and limited streamer tubes interleaved with layers of steel 
and brass. 

The results presented here are obtained using a to- 
tal integrated luminosity of 347.5 fb^^. Monte Carlo 
(MC) simulation samples of T(4S') decays, charm, 
and light quark pairs from continuum, equivalent to 
3.3, 1.7, and 1.1 times the data statistics, respectively, 
and have been generated using Geant4 [38]. These sam- 
ples are used mainly to evaluate background components. 
Quark fragmentation in continuum events is described 
using the JETSET package [39]. The MC distributions 
are rescaled to the data sample luminosity, using the ex- 
pected cross sections of the diff'erent components : 1.3 
nb for cc, 0.525 nb for B+B' and and 2.09 nb 

for light uu, dd, and ss quark events. Dedicated sam- 
ples of pure signal events, equivalent to 4.5 times the 
data statistics, are used to correct measurements for ef- 
ficiency and finite resolution effects. Radiative decays 
{D+ K-n+e+iy^j) are modeled by PHOTOS [40]. 
Events with a decaying into K~n^TT^ are also re- 
constructed in data and simulation. This control sample 
is used to adjust the c-quark fragmentation distribution 
and the kinematic characteristics of particles accompa- 
nying the meson in order to better match the data. 
It is used also to measure the reconstruction accuracy 
of the missing neutrino momentum. Other samples with 
a £>°, a £)*+, or a Df meson exclusively reconstructed 
are used to define corrections on production character- 
istics of charm mesons and accompanying particles that 



contribute to the background. 

VI. ANALYSIS METHOD 

Candidate signal events are isolated from T(45) and 
continuum events using variables combined into two 
Fisher discriminants, tuned to suppress T(4S') and con- 
tinuum background events, respectively. Several differ- 
ences between distributions of quantities entering in the 
analysis, in data and simulation, are measured and cor- 
rected using dedicated event samples. 



A. Signal Selection 

The approach used to reconstruct mesons decaying 
into K^TT+e+z^e is similar to that used in previous analy- 
ses studying D" R-e+i^e [41] and D+ K+R-e+Ue 
[42]. Charged and neutral particles are boosted to the 
CM system and the event thrust axis is determined. A 
plane perpendicular to this axis is used to define two 
hemispheres. 

Signal candidates are extracted from a sample of events 
already enriched in charm semileptonic decays. Criteria 
applied for first enriching selection are: 

• an existence of a positron candidate with a momen- 
tum larger than 0.5 GeV/c in the CM frame, to 
eliminate most of light quark events. Positron can- 
didates are accepted based on a tight identification 
selection with a pion misidentified as an electron or 
a positron below one per mill; 

• a value of R2 > 0.2, R2 being the ratio between 
second- and zeroth-order Fox- Wolfram moments 
[43], to decrease the contribution from B decays; 

• a minimum value for the invariant mass of the par- 
ticles in the event hemisphere opposite to the elec- 
tron candidate, niopp > 0.5 GeV/c^, to reject lepton 
pairs and two-photon events; 

• the invariant mass of the system formed by the 
positron and the most energetic particle in the can- 
didate hemisphere, rritag > 0.13 GeV/c^, to remove 
events where the lepton is the only particle in its 
hemisphere. 

A candidate is a positron, a charged kaon, and a 
charged pion present in the same hemisphere. A vertex 
is formed using these three tracks, and the corresponding 
probability larger than 10"'' are kept. The value of 
this probability is used in the following with other infor- 
mation to reject background events. 

All other tracks in the hemisphere are defined as "spec- 
tators". They most probably originate from the beam 
interaction point and are emitted during hadronization 
of the created c and c quarks. The "leading" particle is 
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the spectator particle having the highest momentum. In- 
formation from the spectator system is used to decrease 
the contribution from the combinatorial background. As 
charm hadrons take a large fraction of the charm quark 
energy, charm decay products have, on average, higher 
energies than spectator particles. 

To estimate the neutrino momentum, the (i4r~7r+e+z^e) 
system is constrained to the D+ mass. In this fit, esti- 
mates of the direction and of the neutrino energy are 
included from measurements obtained from all tracks reg- 
istered in the event. The £>+ direction estimate is taken 
as the direction of the vector opposite to the momentum 
sum of all reconstructed particles but the kaon, the pion, 
and the positron. The neutrino energy is evaluated by 
subtracting from the hemisphere energy the energy of re- 
constructed particles contained in that hemisphere. The 
energy of each hemisphere is evaluated by considciring 
that the total CM energy is distributed between two ob- 
jects of mass corresponding to the measured hemisphere 
masses [44]. As a Z)+ is expected to be present in the 
analyzed hemisphere and as at least a D meson is pro- 
duced in the opposite hemisphere, minimum values for 
hemisphere masses arc imposed. 

For a hemisphere i, with the index of the other hemi- 



sphere noted as j, the energy i^^^^ and the mass m\^^ 
are defined as: 



hem 2 

-J V / ■ 



2,(i) 



2,(j) 



(measured), m£)) 



(33) 



The missing energy in a hemisphere is the difference 
between the hemisphere energy and the sum of the energy 
of the particles contained in this hemisphere [E™^ = 
Ehem — X]r=i" -^i)- -"-^ ^ given collision, some of the re- 
sulting particles might take a path close to the beam line, 
being therefore undetected. In such cases, as one uses all 
reconstructed particles in an event to estimate the D me- 
son direction, this direction is poorly determined. These 
events are removed by only accepting those in which the 
cosine of the angle between the thrust axis and the beam 
line, cos{6thrust) , is smaller than 0.7. In cases where there 
is a loss of a large fraction of the energy contained in the 
opposite hemisphere, the reconstruction of the D is also 
damaged. To minimize the impact of these cases, events 
with a missing energy in the opposite hemisphere greater 
than 3 GeV are rejected. 

The mass-constrained fit also requires estimates of the 
uncertainties on the angles defining the direction and 
on the missing energy must also be provided. These esti- 
mates are parameterized versus the missing energy in the 
opposite hemisphere which is used to quantify the quality 
of the reconstruction in a given event. Parameterizations 
of these uncertainties are obtained in data and in simula- 
tion using events with a reconstructed — >■ K~Tr~^Tr~^ , 
for which we can compare the measured Z)+ direction 
with its estimate using the algorithm employed for the 
analyzed semileptonic decay channel. — > K^tt^tt^ 



events also allow one to control the missing energy es- 
timate and its uncertainty. Corresponding distributions 
obtained in data and with simulated events are given in 
Fig. 4. These distributions are similar, and the remaining 
differences are corrected as explained in Section VIC 2. 

Typical values for the reconstruction accuracy of kine- 
matic variables, obtained by fitting the sum of two Gaus- 
sian distributions for each variable, are given in Ta- 
ble VI. These values are only indicative as the match- 
ing of reconstructed-to-generated kinematic variables of 
events in five dimensions is included, event-by-event, in 
the fitting procedure. 



TABLE VI: Expected resolutions for the five variables. They 
are obtained by fitting the distributions to the sum of two 
Gaussian functions. The fraction of events fitted in the broad 
component is given in the last column. 

variable cri a2 fraction of events 
in broadest Gaussian 



cos^e 0.068 0.325 
cos6'if 0.145 0.5 
X (rad) 0.223 1.174 
(GeV') 0.081 0.264 
TUKn (GeV/c^) 0.0027 0.010 



0.139 

0.135 
0.135 
0.205 
0.032 



B. Background rejection 

Background events arise from T(4S') decays and 
hadronic events from the continuum. Three variables 
are used to decrease the contribution from BB events: 
i?2, the total charged and neutral multiplicity, and the 
sphericity of the system of particles produced in the event 
hemisphere opposite to the candidate. These variables 
use topological differences between events with B decays 
and events with cc fragmentation. The particle distribu- 
tion in T(45) decay events tends to be isotropic as the 
B mesons are heavy and produced near threshold, while 
the distribution in cc events is jet-like as the CM en- 
ergy is well above the charm threshold. These variables 
arc combined linearly in a Fisher discriminant [45], Fht, 
and corresponding distributions are given in Fig. 5. The 
requirement Fbh > retains 70% of signal and 15% of 
Bi3-background events. 

Background events from the continuum arise mainly 
from charm particles, as requiring an electron and a kaon 
reduces the contribution from light-quark flavors to a low 
level. Because charm hadrons take a large fraction of the 
charm quark energy, charm decay products have higher 
average energies and different angular distributions (rel- 
ative to the thrust axis or to the D direction) as com- 
pared to other particles in the hemisphere, emitted from 
the hadronization of the c and c quarks. The me- 
son decays also at a measurable distance from the beam 
interaction point, whereas background event candidates 
contain usually a pion from fragmentation. Therefore, to 
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FIG. 4: (color online) Distributions of the difference (left) 
between reconstructed and expected values, in the CM frame, 
for direction angles [6, (j)) and for the missing energy in 
the candidate hemisphere. These distributions are normalized 
to the same number of entries. The D"*" is reconstructed in the 
K~Tr'^Tr'^ decay channel. Distributions on the right display 
the relative difference between the histograms given on the 
left. 



decrease the amount of background from fragmentation 
particles in cc events, the following variables are used: 

• the spectator system mass; 

• the momentum of the leading spectator track; 

• a quantity derived from the probability of the 

mass-constrained fit; 

• a quantity derived from the vertex fit probability 
of the K , TT and e trajectories; 

• the value of the momentum after the £)+ mass- 
constrained fit; 

• the significance of the flight length of the £>+ from 
the beam interaction point until its decay point; 

• the ratio between the significances of the distance 
of the pion trajectory to the decay position and 
to the beam interaction point. 

Several of these variables are transformed such that dis- 
tributions of resulting (derived) quantities have a bell- 
like shape. These seven variables are combined linearly 
into a Fisher discriminant variable (Fee) and the corre- 
sponding distribution is given in Fig. 6; events are kept 




FIG. 5; (color online) Distributions of Ftt for signal and for 
T(4S') background events. The two distributions are normal- 
ized to the same number of entries. 
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FIG. 6: (color online) Fisher discriminant variable Fee dis- 
tribution for charm background and signal events. The two 
distributions are normalized to the same number of entries. 



for values above 0.5. This selection retains 40% of sig- 
nal events that were kept by the previous selection re- 
quirement on Fhb and rejects 94% of the remaining back- 
ground. About 244 X 10"^ signal events are selected with a 
ratio S/B = 2.3. In the mass region of the X*(892)° this 
ratio increases to 4.6. The average efficiency for signal 
is 2.9% and is uniform when projected onto individual 
kinematic variables. A loss of efficiency, induced mainly 
by the requirement of a minimal energy for the positron, 
is observed for negative values of cos^e and at low q"^. 



C. Simulation tuning 

Several event samples are used to correct differences 
between data and simulation. For the remaining T(45) 
decays, the simulation is compared to data as explained 
in Section VI CI. For e+e" — >■ cc events, corrections 
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to the signal sample are different from those to the 
background sample. For signal, events with a recon- 
structed — )■ K~Tr~^TT~^ in data and MC are used. 
These samples allow us to compare the different distri- 
butions of the quantities entering in the definition of 
the Fbi, and Fee discriminant variables. Measured dif- 
ferences are then corrected, as explained below (Section 
VIC 2). These samples are used also to measure the re- 
construction accuracy on the direction and missing en- 
ergy estimates for £)+ — >• K~TT~^e~^i/e. For background 
events (Section VI C 3), the control of the simulation has 
to be extended to D*+ and Df production and to 
their accompanying charged mesons. Additional sam- 
ples with a reconstructed exclusive decay of the corre- 
sponding charm mesons are used. Corrections are ap- 
plied also on the semileptonic decay models such that 
they agree with recent measurements. Effects of these 
corrections are verified using wrong sign events (Section 
VIC 4), which are used also to correct for the produc- 
tion fractions of charged and neutral Z?-mesons. Finally, 
absolute mass measurement capabilities of the detector 
and the mass resolution are verified (Section VI C 5) us- 
ing K~-K^ and Z)+ — )• i^~7r+7r+ decay channels. 



1. Background from T{AS) decays 

The distribution of a given variable for events from the 
remaining T{AS) BB background is obtained by com- 
paring corresponding distributions for events registered 
at the T(AS) resonance and 40 MeV below. Compared 
with expectations from simulated events in Fig. 7, distri- 
butions versus the kinematic variables agree reasonably 
well in shape, within statistics, but the simulation needs 
to be scaled by 1.7 ± 0.2. A similar effect was measured 
also in a previous analysis of the — )• K~K'^e^Ve de- 
cay channel [42]. 



2. Simulation tuning of signal events 

Events with a reconstructed £>+ — >■ ii'~7r+7r+ candi- 
date are used to correct the simulation of several quanti- 
ties which contribute to the K~iT'^e'^Vf, event reconstruc- 
tion. 

Using the K~'k'^tt^ mass distribution, a signal region, 
between 1.849 and 1.889 GeV/c^, and two sidebands 
([1.798, 1.838] and [1.900, 1.940] GeV/c^), are defined. 
A distribution of a given variable is obtained by subtract- 
ing from the corresponding distribution of events in the 
signal region half the content of those from sidebands. 
This approach is referred to as sideband subtraction in 
the following. It is verified with simulated events that 
distributions obtained in this way agree with those ex- 
pected from true signal events. 

a. control of the c — )• production mechanism: 
the Fisher discriminants and Fee arc functions of sev- 
eral variables, listed in Section VI B, which have distribu- 
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FIG. 7: Ratio (data/MC) distribution for T{AS) decays versus 
each of the five kinematic variables. The dotted line corre- 
sponds to data/MC = 1.7. 



tions that may differ between data and simulation. For 
a given variable, weights arc computed from the ratio of 
normalized distributions measured in data and simula- 
tion. This procedure is repeated, iteratively, considering 
the various variables, until corresponding projected dis- 
tributions are similar to those obtained in data. There 
are remaining differences between data and simulation 
coming from correlations between variables. To mini- 
mize their contribution, the energy spectrum of — >■ 
K'-K^-K^ is weighted in data and simulation to be similar 
to the spectrum of semileptonic signal events. 

We have performed another determination of the cor- 
rections without requiring that these two encirgy spectra 
are similar. Differences between the fitted parameters 
obtained using the two sets of corrections axe taken as 
systematic uncertainties. 

h. control of the direction and missing energy 
measurements: the direction of a fully reconstructed 
D"*" — >■ /ir^7r+7r+ decay is accurately measured and one 
can therefore compare the values of the two angles, defin- 
ing its direction, with those obtained when using all 
particles present in the event except those attributed 
to the decay signal candidate. The latter procedure 
is used to estimate of the D'^ direction for the decay 
Z)+ — > X^TT+e+i^e- Distributions of the difference be- 
tween angles measured with the two methods give the 
corresponding angular resolutions. This event sample al- 
lows also one to compare the missing energy measured 
in the D'^ hemisphere and in the opposite hemisphere 
for data and simulated events. These estimates for the 
D'^ direction and momentum, and their corresponding 
uncertainties are used in a mass-constrained fit. 

For this study, differences between data and simula- 
tion in the c — >■ fragmentation characteristics are 
corrected as explained in the previous paragraph. Global 
cuts similar to those applied for the ^^"Tr+e+i^e 
analysis are used such that the topology of — >• 
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K~TT^Tr~^ selected events is as close as possible to that of 
semileptonic events. Comparisons between angular res- 
olutions measured in data and simulation indicate that 
the ratio data/MC is 1.1 in the tails of the distributions 
(Fig. 4). Corresponding distributions for the missing en- 
ergy measured in the signal hemisphere; {E'^°-^l), in data 
and simulation, show that these distributions have an 
offset of about 100 MeV/c^ (Fig. 4) which corresponds 
to energy escaping detection even in absence of neutri- 
nos. To evaluate the neutrino energy in D+ semileptonic 
decays this bias is corrected on average. 

The difference between the exact and estimated values 
of the two angles and missing energy is measured ver- 
sus the value of the missing energy in the opposite event 
hemisphere {E°^f^^ ). This last quantity provides an es- 
timate of the quality of the energy reconstruction for a 
given event. In each slice of E'^f^^ , a Gaussian distri- 
bution is fitted and corresponding values of the average 
and standard deviation are measured. As expected, the 
resolution gets worse when E'^f^g increases. These val- 
ues are used as estimates for the bias and resolution for 
the considered variable. Fitted uncertainties are slightly 
higher in data than in the simulation. From these mea- 
surements, a correction and a smearing are defined as a 
function of -E^? ^ . They are applied to simulated event 
estimates of 6*, (f) and E^^^ . This additional smearing 
is very small for the D"*" direction determination and is 
typically ~ 100 MeV on the missing energy estimate. 

After applying corrections, the resolution on simulated 
events becomes slightly worse than in data. When eval- 
uating systematic uncertainties we have used the total 
deviation of fitted parameters obtained when applying 
or not applying the corrections. 



3. Simulation tuning of charm background events from 
continuum 

As the main source of background originates from track 
combinations in which particles are from a charm me- 
son decay, and others from hadronization, it is necessary 
to verify that the fragmentation of a charm quark into 
a charm meson and that the production characteristics 
of charged particles accompanying the charm meson are 
similar in data and in simulation. 

In addition, most background events contain a lepton 
from a charm hadron semileptonic decay. The simula- 
tion of these decays is done using the ISGW2 model 
[46], which does not agree with recent measurements [41], 
therefore all simulated decay distributions are corrected. 

a. Corrections on charm quark hadronization: for 
this purpose, distributions obtained in data and MC are 
compared. We study the event shape variables that en- 
ter in the Fisher discriminant F},i, and for variables en- 
tering into Fee, apart from probability of the mass- 
constrained fit which is peculiar to the analyzed 
semileptonic decay channel. Production characteristics 
of charged pions and kaons emitted during the charm 



quark fragmentation, are also measured, and their rate, 
momentum, and angle distribution relative to the sim- 
ulated D direction are corrected. These corrections are 
obtained separately for particles having the same or the 
opposite charge relative to the charm quark forming the 
D hadron. Corrections consist of a weight applied to each 
simulated event. This weight is obtained iteratively, cor- 
recting in turn each of the considered distributions. Mea- 
surements are done for (vetoing from _D*+ 
decays) and for £)+. For £)+ mesons, only the corre- 
sponding c-quark fragmentation distribution is corrected. 

h. Correction of D semileptonic decay form factors: 
by default, D semileptonic decays are generated in Evt- 
Gen [47] using the ISGW2 decay model which does not 
reproduce present measurements (this was shown for in- 
stance in the BaBAR analysis of D*+ -J> D'^Tr+,D° -> 
K~e'^Uf. [41]). Events are weighted such that they cor- 
respond to hadronic form factors behaving according to 
the single pole parameterization as in Eq. (24). 

For decay processes of the type D Pevg, where 
P is a pseudoscalar meson, the weight is proportional 
to the square of the ratio between the correspond- 
ing hadronic form factors, and the total decay branch- 
ing fraction remains unchanged after the transforma- 
tion. For all Cabibbo-favored decays a pole mass value 
equal to 1.893 GeV/c^ [41] is used whereas for Cabibbo- 
suppressed decays 1.9 GeV/c^ [48] is taken. This value 
of the pole mass is used also for Dg semileptonic de- 
cays into a pseudoscalar meson. For decay processes 
of the type D Vei^e, (V P1P2), where P and V 
are respectively pseudoscalar and vector mesons, correc- 
tions depend on the mass of the hadronic system, and 
on q^, cosOe, cos 9k and x- They are evaluated it- 
eratively using projections of the differential decay rate 
versus these variables, as obtained in EvtGen and in a 
simulation which contains the expected distribution. To 
account for correlations between these variables, once dis- 
tributions agree in projection, binned distributions over 
the five dimensional space are compared and a weight 
is measured in each bin. For Cabibbo-allowed decays, 
events are distributed over 2800 bins, similar to those 
defined in Section VI D; 243 bins are used for Cabibbo- 
suppressed decays. Apart for the resonance mass and 
width which are different for each decay channel, the 
same values, given in Table VII, are used for the other 
parameters which determine the differential decay rate. 

For decay channels D — >■ Kwe'^Ve an S-wave compo- 
nent is added with the same characteristics as in the 
present measurements. Other decay channels included 
in EvtGen [47] and contributing to this same final state, 
such as a constant amplitude and the _ft'2(1430)*^ compo- 
nents, are removed as they are not observed in data. 

All branching fractions used in the simulation agree 
within uncertainties with the current measurements [16] 
(apart for D — ^ nC^Vg, which is then rescaled). Only 
the shapes of charm semileptonic decay distributions are 
corrected. 

Systematic uncertainties related to these corrections 
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FIG. 8: (color online) Distributions of the five dynamical variables for wrong-sign events in data (black dots) and MC (his- 
tograms), after all corrections. From top to bottom the background components displayed in the stacked histograms are: 
cc [Ds, , D'^), uds, and BB events respectively. In the lower row, distributions of the ratio data/MC for upper row plots 
are given. 



are estimated by varying separately each parameter ac- 
cording to its expected uncertainty, given in Table VII. 



4- Wrong sign event analysis 



TABLE VII: Central values and variation range for the var- 
ious parameters which determine the differential decay rate 
in D — > P/Ve^Ve decays, used to correct the simulation and 
to evaluate corresponding systematic uncertainties. The form 
factors Ai(q^), A2{q'^) and V(q^) and the mass parameters 
niA and mv are defined in Eq. (24). 



parameter 



central 
value 



variation 
interval 



mpole{D°'' 



ire^Ue) 



1.893 GeV/c^ 
1.9 GeV/c^ 



r)lri'e+Ve) 1.9 GeV/c 



±30 MeV/c^ 
±100 MeV/c^ 
±100 MeV/c^ 



r2 = |^2(0)|/|Ai(0)| 

rv = \vm/\Am\ 

ruA 
mv 
rsw 



0.80 
1.50 

2.5 GeV/c^ 



±0.05 
±0.05 
±0.3 GeV/c^ 



2.1 GeV/c' ±0.2 GeV/c" 
3.0 (GeV/c)-i ±0.3 (GeV/c~') 



Wrong-sign (WS) events of the type K~TT~e^ are used 
to verify if corrections applied to the simulation improve 
the agreement with data, because the origin of these 
events is quite similar to that of the background con- 
tributing in right-sign i^~7r+e+ (RS) events. The ra- 
tio between the measured and expected number of WS 
events is 0.950 ±0.005. In RS events the number of back- 
ground candidates is a free parameter in the fit. 



At this point corrections have been evaluated sepa- 
rately for charged and neutral D mesons. As the two 
charged states correspond to background distributions 
having different shapes, it is also possible to correct for 
their relative contributions. We improve the agreement 
with data by increasing the fraction of events with a 
meson in MC by 4% and correspondingly decreasing the 
fraction of D'^ by 5%. After corrections, projected distri- 
butions of the five kinematic variables obtained in data 
and simulation are given in Fig. 8. 
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5. Absolute mass scale. 



The absolute mass measurement is verified using exclu- 
sive reconstruction of charm mesons in data and simula- 
tion. For candidate events D*+ £'°7r+, £)° K~Tr~^, 
the mean and RMS values of the mass distribu- 
tion arc measured from a fit of the sum to a Gaus- 
sian distribution for the signal and a first order polyno- 
mial for the background. The mass reconstructed 
in simulation is very close to expectation, A^'-^ = 
(-0.07 ± 0.01) MeV/c^ whereas in data it differs by 
A^*" = (-1.07 ±0.17) MeV/c2. Here A„,, is the dif- 
ference between the reconstructed and the exact or the 
world average mass values when analyzing MC or data 
respectively. The uncertainty quoted for Af"'" is from 
Ref. [16]. To correct for this effect the momentum (p) of 
each track in data, measured in the laboratory frame, is 
increased by an amount: A^"*" = 0.7 x lQ~^p. The stan- 
dard deviation of the Gaussian fitted on the D° signal 
is slightly smaller in simulation, (7.25 ± 0.01) MeV/c^, 
than in data, (7.39 ± 0.01) McV/c^. The difference be- 
tween the widths of reconstructed D'^ signals in the two 
samples, is measured versus the transverse momentum of 
the tracks emitted in the decay. In simulation, the mea- 
sured transverse momenta of the tracks are smeared to 
correct for this difference. 

Having applied these corrections, D+ mass distribu- 
tions, for the decay D+ ii'~7r+7r+ obtained in data 
and simulation are comparcxl. The standard deviation of 
the fitted Gaussian distribution on signal is now similar 
in data and simulation. The reconstructed £>+ mass is 
higher by 0.23 MeV/c^ in simulation (on which no cor- 
rection was applied) and by 0.32 MeV/c^ in data. These 
remaining differences are not corrected and included as 
uncertainties. 



D. Fitting procedure 

A binned distribution of data events is analyzed. The 
expected number of events in each bin depends on signal 
and background estimates and the former is a function 
of the values of the fitted parameters. 

We perform a minimization of a negative log-likelihood 
distribution. This distribution has two parts. One cor- 
responds to the comparison between measured and ex- 
pected number of evcmts in bins which span the five di- 
mensional space of the differential decay rate. The other 
part uses the distribution of the values of the Fisher dis- 
criminant variable Fee to measure the fraction of back- 
ground events. 

There are respectively 5, 5 and 4 equal size bins for the 
variables x, cos 9k and cos Be- For and rnxw we use 
respectively 4 and 7 bins of different size such that they 
contain approximately the same number of signal events. 
There are 2800 bins (iVbins) in total. 

The likelihood expression is: 



Whins 




n ^("datal'^Mc) 


i=0 




A^data 

n 







Nhkg 



X Pdfi 



bkg 



(34) 



where n^ata number of data events in bin i and nj^^,-, 
is the sum of MC estimates for signal and background 
events in the same bin. -P('"-datal"'Mc) ^he Poisson 
probability for having n^^^^g^ events in bin i where 71^^! 
events are expected, on average, where: 



+ 



H^tot 



(Ao,A) 



^Din z 
events 

E 

A^bkg 

bkff 



all bins 
Bvents 

E 

j=0 



Slg 



W^^f(Ao,A)W^,(Ao 



N. 



bkg' 



(35) 



Wj{Xo) 



The summation to determine Ti^c extends over all gen- 
erated signal events which are reconstructed in bin i. The 
terms A and Aq are, respectively, the values of parame- 
ters used in the fit and those used to produce simulated 
events. Wj (A) is the value of the expression for the decay 
rate (see Eq. (16)) for event j using the set of parame- 
ters (A). In these expressions, generated values of the 
kinematic variables are used. Cj is the weight applied to 
each signal event to correct for differences between data 
and simulation. It is left unchanged during the fit. W^^.^ 
is the estimated number of background events in bin i 
given by the simulation, corrected for measured differ- 
ences with data, as explained in Section VIC. W^kg 
the estimated total number of background events. 

A^sig and A^bkg are respectively the total number of 
signal and background events fitted in the data sample 
which contains iVdata events, pdf^^^ and pdf^^ are the 
probability density functions for signal and badcground, 
respectively, evaluated at the value of the Fee variable for 
event j. The following expressions are used : 



pdfsig{Fec) = Kig \ C2 X exp 



-{Fee-Cof 



C5 X exp 



pdfhkgiFee) = M,kg \ exp 
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-{Fee-Csr 

2cl 

■ 4 

Y^di^Feef 
.i=0 



(36) 



and values of the corresponding parameters co-5 and 
(io-4 are determined from fits to binned distributions of 
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Fee in simulated signal and background samples. AQig 
and A/'bkg are normalization factors. In Fig. 9 these two 
distributions are drawn to illustrate their different be- 
havior versus the values of Fee for signal and background 
events. As expected, the pdfhkg distribution has higher 
values at low Fee than the corresponding distribution for 
signal. 



1. Background smoothing 

As the statistics of simulated background events for 
the charm continuum is only 1.6 times the data, biases 
appear in the determination of the fit parameters if we 
use simply, as estimates for background in each bin, the 
actual values obtained from the MC. Using a parameter- 
ized event generator, this effect is measured using distri- 
butions of the difference between the fitted and exact val- 
ues of a parameter divided by its fitted uncertainty (pull 
distributions). To reduce these biases, a smoothing [49] 
of the background distribution is performed. It consists 
of distributing the contribution of each event, in each di- 
mension, according to a Gaussian distribution. In this 
procedure correlations between variables are neglected. 
To account for boundary effects, the dataset is refiected 
about each boundary, x is essentially uncorrelated with 
all other variables and in particular with cos 9i. There- 
fore, for each bin in (m, and cos 9k), a smoothing of 
the X and cos 0i distributions is done in the hypothesis 
that these two variables are independent. 



0.15 




FIG. 9: (color online) Probability density functions for signal 
(red dashed line) and background (blue full line) events versus 
the values of the discriminant variable Fee- 



VII. D+ K*°e+u, HADRONIC FORM FACTOR 
MEASUREMENTS 

We first consider a signal made of the if* (892)" and S- 
wave components. Using the LASS parameterization of 
the S'-wave phase versus the Ktt mass (Eq. (10)), values 



of the following quantities (quoted in Table VIII second 
column) are obtained from a fit to data : 

• parameters of the if* (892)° Breit-Wi gner distri- 
bution: 771^^.(892), r/f*(892)' ^^'^ "i^BW, the Blatt- 
Weisskopf parameter; 

• parameters of the K e^Ve hadronic form 
factors: r2, ry, and m^. The parameter my which 
determines the variation of the vector form fac- 
tor is fixed to 2.0 GeV/c^; 

• parameters which define the 5- wave component: rg 
and r'^p for the amplitude (Eq. (28)), a]j%Q and 

1 /2 

BG t^*^ phase (Eq. (8)); 

• and finally the total numbers of signal and back- 
ground events, N^ig and N^kg- 

1 /2 

Apart from the effective range parameter, hg ^q, all 
other quantities are accurately measured. Values for 
the S-wave parameters depend on the parameterization 
used for the P-wave and as the LASS experiment in- 
cludes a if* (1410) and other components one cannot 
directly compare our results on a^gQ and b^J^u with 
those of LASS. We have obtained the first measurement 
for niA which gives the q^ variation of the axial vector 
hadronic form factors. Using the values of fitted parame- 
ters and integrating the corresponding differential decay 
rates, fractions of S- and P-wave are given in the second 
column of Table IX. 

Projected distributions, versus the five variables, ob- 
tained in data and from the S-wave + K (892)° fit resuh 
are displayed in Fig. 10. The total x^ of this fit is 2914 for 
2787 degrees of freedom which corresponds to the prob- 
ability of 4.6%. Fit results including the if* (1410)° and 
D-wave are discussed in Section VIII. 

A. Systematic uncertainties 

The systematic uncertainty on each fitted parameter 
(x) is defined as the difference between the fit results 
in nominal conditions (a;[0]) and those obtained, (^^[j]), 
after changing a variable or a condition (i) by an amount 
which corresponds to an estimate of the uncertainty in 
the determination of this quantity: 

Ax ^ x[0] - x\i]. (37) 

Values are given in Table X. Some of the corrections 
induce a variation on the Fee distributions for signal or 
background which are therefore reevaluated. 

1. Signal production and decay 

a. Corrections of distributions of Fisher input vari- 
ables (I): the signal control sample is corrected as ex- 
plained in Section VIC 2. The corresponding systematic 
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TABLE VIII: Values of fitted parameters assuming that the final state consists of a sum of 5- wave and K (892)" components 
(second column), includes the K (1410)" in the P-wavc (third column) and a _D-wave (last column). The variation of the 
S-wave phase versus the K-k mass is parameterized according to Eq. (10) whereas the 5*- wave amplitude is parameterized as in 
Eq. (28). Fit results including the K (1410)° are discussed in Section VIII A. Values given in the third column of this Table 
are the central results of this analysis. 



variable 


S + K*{892)° 


5 ± if* (892)" 


5 ± if* (892)° 
K (1410) + D 






K (1410) 


mx.(892)(McV/c2) 
r?c*(892)(MeV/c2) 


894.77 ± 0.08 


895.43 ± 0.21 


895.27 ±0.15 


45.78 ± 0.23 


46.48 ±0.31 


46.38 ± 0.26 


rsM'(GeV/c)-i 

mA(GcV/c2) 


3.71 ±0.22 


2.13 ±0.48 


2.31 ± 0.20 


2.65 ±0.10 


2.63 ±0.10 


2.58 ±0.09 


TV 


1.458 ±0.016 


1.463 ±0.017 


1.471 ±0.016 




0.804 ± 0.020 


0.801 ±0.020 


0.786 ± 0.020 


rs(GeV)-i 


-0.470 ± 0.032 


-0.497 ±0.029 


-0.548 ± 0.027 




0.17 ±0.08 


0.14 ±0.06 


0.03 ± 0.06 


4(L(GcV/c)-i 


1.82 ± 0.14 


2.18 ± 0.14 


2.10 ±0.10 




-1.66 ±0.65 


1.76 fixed 


1.76 fixed 


^if* (1410)0 




0.074 ± 0.016 


0.052 ± 0.013 


<5if- (1410)0 (degree) 
ri3(GeV)-* 




8.3 ± 13.0 


fixed 






0.78 ±0.18 


5d (degree) 






fixed 


■A^ sig 


243850 ± 699 


243219 ± 713 


243521 ± 688 


Nbkg 


107370 ± 593 


108001 ± 613 


107699 ± 583 


Fit probability 


4.6% 


6.4% 


8.8% 



TABLE IX: Fractions for signal components assuming that the final state consists of a sum of S-wave and K (892)" components 
(second column), including the K (1410)" in the P-wave (third column) and a D-wave (last column). In the second and third 
cases, the sum of the fractions for the two K does not correspond exactly to the total P-wave fraction because of interference. 



Component 


S + K* (892)° 


S'±i?*(892)" 


S + i?*(892)" 




(%) 


±:^*(1410)"(%) 


i?*(1410)" ±D(%) 


S-wave 


5.62 ±0.14 ±0.13 


5.79 ±0.16 ±0.15 


5.69 ±0.16 ±0.15 


P-wave 


94.38 


94.21 


94.12 


F*(892)" 


94.38 


94.11 ±0.74 ±0.75 


94.41 ±0.15 ±0.20 


Z*(1410)" 





0.33 ±0.13 ±0.19 


0.16 ±0.08 ±0.14 


D-wave 








0.19 ±0.09 ±0.09 



uncertainty is obtained by defining new event weights 
without taking into account that the momentum distri- 
bution of reconstructed D mesons is different in hadronic 
and in semileptonic samples. 



c. Particle identification efficiencies (III): the sys- 
tematic uncertainty is estimated by not correcting for 
remaining differences between data and MC on particle 
identification. 



b. Simulation of radiative events (II): most of ra- 
diative events correspond to radiation from the charged 
lepton, although a non- negligible fraction comes from ra- 
diation of the K*(892)° decay products. In D ^ Peue, 
by comparing two generators (PHOTOS [40] and KLOR 
[50]), the CLEO-c collaboration has used a variation of 
16% to evaluate corresponding systematic uncertainties 
[51]. We have increased the fraction of radiative events 
(simulated by PHOTOS) by 30% (keeping constant the 
total number of events) and obtained the corresponding 
variations on fitted parameters. 



d. Estimates of the values and uncertainties for the 
D direction and missing energy (IV): in Section VIC 2 
it is observed that estimates of the £>+ direction and 
energy arc more accurate in the simulation than in data. 
After applying smearing corrections, the result of this 
comparison is reversed. The corresponding systematic 
uncertainty is equal to the difference on fitted parameters 
obtained with and without smearing. 
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FIG. 10: (color online) Projections of data (black dots) and of the fitted MC distribution (histograms) versus each of the five 
kinematic variables. The signal contains an S-wave and a K (892)° components. From top to bottom the fitted background 
components displayed in the stacked histograms are: cc, uds, and BB events respectively. In the lower row, distributions of 
the ratio data/MC for upper row plots are given. 



2. BB background correction (V) 

The number of remaining BB background events ex- 
pected from simulation is rescaled by 1.7 ± 0.2 (see Sec- 
tion VIC 1). The uncertainty on this quantity is used to 
evaluate corresponding systematic uncertainties. 



3. Corrections to the cc background 

a. Fragmentation associated systematic uncertainties 
(VI): after applying corrections explained in Section 
VI C 3 a, the remaining differences between data and sim- 
ulation for the considered distributions are five times 
smaller. Therefore, 20% of the full difference measured 
before applying corrections is used as the systematic un- 
certainty. 

b. Form-factor correction systematics (VII): corre- 
sponding systematic uncertainties depend on uncertain- 
ties on parameters used to model the differential semilep- 
tonic decay rate of the various charm mesons (see Section 
VIC 3 b). 



c. Hadronization- associated systematic uncertainties 
(VIII): using WS events, it is found in Section VIC 4 
that the agreement between data and simulation im- 
proves by changing the hadronization fraction of the dif- 
ferent charm mesons. Corresponding variations of rela- 
tive hadronization fractions are compatible with current 
experimental uncertainties on these quantities. The cor- 
responding systematic uncertainty is obtained by not ap- 
plying these corrections. 

4-. Fitting procedure 

a. Background Smoothing (IX): the MC back- 
ground distribution is smoothed as explained in Section 
VI D 1. The evaluation of the associated systematic un- 
certainty is performed by measuring with simulations 
based on parameterized distributions, the dispersion of 
displacements of the fitted quantities when the smooth- 
ing is or is not applied in a given experiment. It is verified 
that uncertainties on the values of the two parameters 
used in the smoothing have negligible contributions to 
the resulting uncertainty. 
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TABLE X: Systematic uncertainties on parameters fitted using the S-wave and -f!r*(892)'' model, expressed as {x[0]—x[i])/(7stat' 
(I) uncertainty associated with the tuning of the signal control sample, (II) fraction of radiative signal events increased by 30%, 
(III) no PID corrections on electron or kaon in MC signal events, (IV) no smearing applied on 6d, 4>d and Emiss for simulated 
signal events, (V) BB background rate lowered by the statistical uncertainty of its determination, (VI) uncertainty associated 
with the tuning of fragmentation in charm background events, (VII) remaining uncertainty on semileptonic decay models in 
charm background events, (VIII) uncertainty associated with c-meson relative fractions, (IX) uncertainty remaining from the 
smoothing of the background distribution, (X) effects from limited statistics in simulation, (XI) variation of parameters that 
were kept constant in the fit and, (XII) absolute mass scale uncertainties. 

variation AMjf.(s92)0 ^^K'{sn2)" Arsiv AniA Ary Ar2 Ars Ar^g^ Aa^g^^ Ab^gf^ ANs ANb 
signal 

I ^OlS 

II -0.36 

III 0.21 

IV 0.29 



BB bkg. 


V 


-0.06 


0.32 


0.09 


0.22 


-0.13 


0.03 


0.30 


0.31 


0.14 


0.30 


-0.09 


0.11 


cc bkg. 


VI 


-0.04 


0.21 


-0.61 


0.10 


-0.08 


0.07 


0.33 


0.32 


0.13 


0.27 


0.06 


-0.08 


VII 


0.53 


0.19 


0.14 


0.16 


0.13 


0.07 


0.10 


0.10 


0.17 


0.19 


0.16 


0.22 


VIII 


0.24 


0.36 


0.11 


-0.49 


0.85 


0.04 


-0.76 


-0.68 


-0.77 


1.02 


0.76 


-0.91 


Fitting procedure 


IX 


0.13 


0.17 


0.25 


0.29 


0.30 


0.25 


0.25 


0.25 


0.32 


0.32 


0.13 


0.13 


X 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


XI 


0.00 


0.00 


0.07 


0.07 


1.15 


0.08 


0.05 


0.05 


1.43 


0.46 


0.01 


0.01 


XII 


-0.93 


-0.06 


0.09 


0.09 


-0.05 


0.04 


0.03 


0.02 


0.07 


-0.05 


0.00 


0.00 


(Tsyst 


1.41 


1.00 


1.06 


1.21 


1.87 


0.97 


1.27 


1.23 


1.87 


1.47 


1.29 


1.48 



-0.16 -0.10 -0.18 0.28 0.18 -0.40 

0.07 0.02 -0.11 0.34 0.10 0.26 

0.13 0.27 0.69 0.78 0.51 0.17 

0.36 0.20 -0.18 0.07 -0.25 0.15 



-0.43 0.02 0.00 -0.36 0.44 

0.20 0.17 0.21 -0.21 0.26 

0.16 0.29 0.17 0.18 0.22 

0.19 -0.31 -0.23 0.57 -0.70 



b. Limited statistics of simulated events (X): fluc- 
tuations of the number of MC events in each bin are not 
included in the likelihood expression, therefore one quan- 
tifies this effect using fits on distributions obtained with a 
parameterized event generator. Pull distributions of fit- 
ted parameters, obtained in similar conditions as in data, 
have an RMS of 1.2. This increase is attributed to the 
limited MC statistics used for the signal (4.5 times the 
data) and, also, from the available statistics used to eval- 
uate the background from e+e~ — > cc continuum events. 
We have included this effect as a systematic uncertainty 
corresponding to 0.7 times the quoted statistical uncer- 
tainty of the fit. It corresponds to the additional fluctu- 
ation needed to obtain a standard deviation of 1.2 of the 
pull distributions. 



6. Absolute mass scale (XII). 

When corrections defined in Section VI C 5 are applied, 
in data and simulation, for the — )• i>sr~7r+e+z^e decay 
channel, the fitted K*(892)'^ mass in data increases by 
0.26 MeV/c^ and its width decreases by 0.12 MeV/c^. 
The uncertainty on the absolute mass measurement of 
the if* (892)° is obtained by noting that a mass variation, 
^data^ of the D reference signal is reduced by a factor of 
four in the K* mass region; this gives: 

a(mK.(892)o) = V0.17^ + 0.23^ ^|;^^^^J^ 



~ 0.07 MeV/c^ 



(38) 



5. Parameters kept constant in the fit (XI). 



In this expression, 0.17 MeV/c^ is the uncertainty on the 
mass [16] and 0.23 MeV/c^ is the difference between 
the reconstructed and exact values of the D+ mass in 

simulation (s(h; Section VIC 5). Uncertainty on the; 
K* width measurement from track resolution effects is 
negligible. 



The signal model has three fixed parameters, the 
vector pole mass my and the mass and width of the 
-ftro(1430) resonance. Corresponding systematic uncer- 
tainties are obtained by varying the values of these pa- 
rameters. For mv a ±100 MeV/c^ variation is used, 
whereas for the other two quantities we take respectively 
±50 MeV/c2 and ±80 MeV/c^ [16]. 



7. Comments on systematic uncertainties 

The total systematic uncertainty is obtained by sum- 
ming in quadrature the various contributions. The main 

systematic uncertainty on ry comes from the assumed 
variation for the parameter my because these two pa- 
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TABLE XI: Systematic uncertainties on parameters fitted using a model for the signal which contains S-w&ve, K* (892)" and 
K (1410)" components, expressed as (x-[0] — x[i\) / a stat- (I) uncertainty associated with the tuning of the signal control sample, 
(II) fraction of radiative signal events increased by 30%, (III) no PID corrections on electron or kaon in MC signal events, 
(IV) no smearing applied on 9d, 4>d and Emiss for sirrmlated signal events, (V) BB background rate lowered by the statistical 
uncertainty of its determination, (VI) uncertainty associated with the tuning of fragmentation in charm background events, 
(VII) remaining uncertainty on semileptonic decay models for background events, (VIII) uncertainty associated with c-meson 
relative fractions, (IX) uncertainty remaining from the smoothing of the background distribution, (X) effects from limited 
statistics in simulation, (XI) variation of parameters that were kept constant in the fit, (XII) uncertainties on absolute mass 
scale. 



variation 


^^K* (892)0 


^'^K* (892)0 


Arsw 


AruA 


Arv 


Ar2 


Ars 


a4^) 




ft:* (1410)0 


^''fi:* (1410)0 


ANs 


AATb 


signal 


(I) 


0.17 


0.05 


-0.23 


-0.22 


-0.31 


0.18 


0.14 


-0.14 


-0.13 


0.23 


-0.19 


-0.39 


0.45 


(11) 


-0.18 


0.06 


-0.01 


-0.14 


-0.36 


0.09 


-0.10 


0.08 


0.05 


-0.08 


-0.08 


-0.23 


0.26 


(III) 


0.02 


0.10 


0.06 


0.70 


0.73 


0.53 


0.14 


0.07 


0.41 


0.08 


0.22 


0.10 


0.12 


(IV) 


-0.13 


0.03 


0.29 


-0.18 


-0.04 


-0.27 


-0.02 


0.02 


-0.17 


-0.18 


0.32 


0.61 


-0.70 


BB bkg. 


(V) 


-0.41 


-0.04 


0.34 


0.26 


0.16 


0.05 


-0.12 


0.12 


0.08 


-0.46 


0.22 


-0.01 


0.02 


cc bkg. 


(VI) 


-0.14 


0.07 


-0.08 


0.13 


0.09 


0.09 


-0.16 


0.14 


-0.01 


-0.24 


-0.03 


0.08 


-0.09 


(VII) 


0.09 


0.08 


0.14 


0.19 


0.14 


0.08 


0.18 


0.15 


0.06 


0.10 


0.11 


0.14 


0.18 


(VIII) 


-0.44 


-0.19 


0.59 


-0.48 


-0.75 


0.04 


0.98 


-0.94 


0.28 


-0.42 


-0.21 


1.03 


-1.23 


Fitting proccHlur(> 


(IX) 


0.13 


0.17 


0.25 


0.29 


0.30 


0.25 


0.25 


0.25 


0.32 


0.30 


0.30 


0.13 


0.13 


(X) 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


(XI) 


0.27 


0.12 


0.29 


0.07 


1.15 


0.08 


0.57 


0.55 


3.25 


0.89 


0.40 


0.09 


0.10 


(XII) 


-0.33 


-0.05 


0.03 


0.09 


0.05 


0.04 


-0.05 


0.03 


0.06 


-0.02 


-0.01 


-0.02 


0.02 




1.08 


0.78 


1.13 


1.24 


1.81 


0.99 


1.40 


1.35 


3.39 


1.39 


1.02 


1.48 


1.69 



ramctcrs arc correlated. Values of the parameters cig 

and b^gQ depend on the mass and width of the Xq(1430) 
because the measured 5-wave phase is the sum of two 
components: a background term and the ArQ(1430). 



VIII. INCLUDING OTHER COMPONENTS 

A contribution to the P-wave from the K (1410)° ra- 
dial excitation was measured by LASS [5] in Kp inter- 
actions at small transfer and in r decays [12]. As is 
discussed in the following, even if the statistical signif- 
icance of a signal at high mass does not reach the level 
to claim an observation, data favor such a contribution 
and a signal containing the X*(892)°, the X*(1410)" and 
an S-wave components is considered as our nominal fit 
to data. 

To compare present results for the S-wave with 
LASS measurements a possible contribution from the 
K (1410)" is included in the signal model. It is parame- 
terized using a similar Breit-Wigner expression as for the 
K (892)° resonance. The L=l form factor components 
are in this case written as: 



J"ii oc {BW + rK-'{i4io)oe'^'<'(^"°>" BW')2V2qHo (39) 
T21 oc (BW^-hr^.(i4io)oe*'^^*("i'')"W)2g(i?+-f i?_) 
J"3i oc {BW + rK-'{i4io)oe'^'''<^'"°>° BW')2q{H+ - H_) 

where BW stands for the IT (892)° Breit-Wi gner distri- 
bution (Eq. (25)) and BW for that of the ;^*(1410)°. 
As the phase space region where this last component 
contributes is scarcely populated (high Ktt mass), this 
analysis is not highly sensitive to the exact shape of 
the resonance. Therefore the Breit-Wigner parameters 
of the K (1410)° (given in Table I) are fixed and only 
the relative strength {rK»(i4io)") ^^id phase ((5x* (1410)0) 

1/2 

are fitted. For the same reason, the value of bg = 

1.76 GeV"^ is fixed to the LASS result (given in Ta- 
ble III). 

A. Results with a j£r*(1410)° contribution included 

Results are presented in Table VIII (third column) us- 
ing the same S'-wavc parameterization as in Section VII. 
They correspond to the central results of this analysis. 

The total value is 2901 and the number of degrees 
of freedom is 2786. This corresponds to a probability of 
6.4%. Systematic uncertainties, evaluated as in Section 
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FIG. 11: (color online) Projections of data (black dots) and of the fitted MC distribution (histograms) versus each of the five 
kinematic variables. The signal contains S'-wave, K (892)° and K (1410)" components. From top to bottom the background 
components displayed in the stacked histograms are: cc, uda, and BB events respectively. In the lower row, distributions of 
the ratio data/MC for upper row plots are given. 
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FIG. 12: (color online) Comparison between measured and fit- 
ted distributions of the values of the Fee discriminant variable. 
Points with error bars correspond to data. The histogram is 
the fitted distribution. It is the sum of a background (blue, 
filled histogram) and signal (hatched) components. 
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FIG. 13: (color online) Comparison between measured and 
fitted K-K mass distributions in the K (892)" region. Results 
of a fit in which the width of the K (892)" meson is fixed 
to 50.3 MeV/c^ (value quoted in 2008 by the Particle Data 
Group) are also given. 



VII A, are given in Table XI. The statistical error matrix 
of fitted parameters, a table showing individual contri- 



bution of sources of systematic uncertainties, which were 
grouped in the entries of Tab. XI labelled III, VII and 
XI, and the full error matrix of systematic uncertainties 
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FIG. 14: (color online) Background subtracted data mass distribution (black full dots) and fit result (red open crosses) for 
the S'-wave and K*{892)° model (left) and the S-wave, K* {892)" and ii:*(1410)° model (right), in the high mass region. Error 
bars correspond to statistical uncertainties only. 



are given in A. Projected distributions versus the five 
variables obtained in data and from the fit result are dis- 
played in Fig. 11. Measured and fitted distributions of 
the values of the F^c discriminant variable are compared 
in Fig. 12. 

The comparison between measured and fitted, back- 
ground subtracted, mass distributions is given in Fig. 13. 
Results of a fit in which the width of the if * (892)° reso- 
nance IS fixed to 50.3 MeV/c^ (the value quoted in 2008 
by the Particle Data Group) are also given. 

Background subtracted projected distributions versus 
niKTT for values higher than 1 GeV/c^, obtained in data 
and using the fit results with and without including the 
^*(1410)", are displayed in Fig. 14. 

The measured fraction of the K (1410)" is compatible 
with the value obtained in r decays [12]. The relative 
phase between the X*(892)° and !K*(1410)" is compati- 
ble with zero, as expected. Values of the hadronic form 

I — *o I 

factor parameters for the decay — > K e^Ve are al- 
most identical with those obtained without including the 
K (1410)". The fitted value for a^J^Q is compatible with 
the result from LASS reported in Table III. 

The total fraction of the S'-wave is compatible with the 
previous value. Fractions for each component are given 
in the third column of Table IX. 

Considering several mass intervals, background sub- 
tracted projected distributions versus the four other vari- 
ables, obtained in data and from the fit results, are dis- 
played in Fig. 15 to 18. 



B. Fit of the A'*(1410)° contribution and of the 
5-wave amplitude and phase 

Fixing the parameters which determine the if* (892)" 
contribution to the values obtained in the previous fit, 
we measure the S'-wave parameters entering in Eq. (28) 
in which the S'-wave phase is assumed to be a constant 
within each of the considered Kt: mass intervals. Values 
of rriKjr which correspond to the center and to half the 
width of each mass interval are given in Table XIII. The 
two parameters which define the K (1410)" are also fit- 
ted. Numbers of signal and background events are fixed 
to their previously determined values. Values of fitted 
parameters are given in Table XIV. 

The variation of the S'-wave phase is given in Fig. 19 
and compared with LASS results and with the result 
found in Section VIII where the S'-wave phase variation 
was parameterized versus the Ktt mass. Systematic un- 
certainties are given in Table XII. 

In Fig. 20 measured values of the S'-wave phase ob- 
tained by various experiments in the elastic region are 
compared. Fig. 20-a is a zoom of Fig. 19. Fig. 20-b to 
-d compare present measurements with those obtained 
in Dalitz plot analyses of the decay — ^ K~7r~^Tr~^. 
For the latter, the S-wave phase is obtained by reference 
to the phase of the amplitude of one of the contribut- 
ing channels in this decay. To draw the different figures 
it is assumed that the phase of the S-wave is equal at 
iTiKTr — 0.67 GeV/c^ to the value given by the fitted pa- 
rameterization on LASS data. It is difficult to draw clear 
conclusions from these comparisons as Dalitz plot anal- 
yses do not provide usually the phase of the I — 1/2 
amplitude alone but the phase for the total S'-wave am- 
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FIG. 15: (color online) Projections of background subtracted data (black dots) and fitted MC signal distributions (hashed 
histogram) versus the four kinematic variables in the mass region between threshold and 800 MeV/c'^. Error bars correspond 
to statistical uncertainties only. The signal contains 5'-wave, K (892)° and K (1410)" components. Lower plots are the ratio 
between data and the fitted signal. 

TABLE Xll: Systematic uncertainties on parameters fitted using a model for the signal which contains S'-wave, K (892)" and 
K (1410)° components in which the K (892)" parameters are fixed, expressed as (x[0] — x[i]) / a stat'-- (1) uncertainty associated 
with the tuning of the signal control sample, (II) fraction of radiative signal events increased by 30%, (III) no PID correction 
on electron or kaon in MC signal events, (IV) no smearing applied on 6d, 4>d and Emiss for simulated signal events, (V) BB 
background rate lowered by the statistical uncertainty of its determination, (VI) uncertainty associated with the tuning of 
fragmentation in charm background events, (VII) remaining uncertainty on semileptonic decay models for background events, 
(Vlll) uncertainty associated with c-meson relative fractions, (IX) uncertainty remaining from the smoothing of the background 
distribution, (X) effects from limited statistics in simulation, (XI) variation of parameters that were kept constant in the fit, 
(Xll) uncertainties on absolute mass scale. 
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cc bkg. 


(VI) 
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0.15 
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1.29 


Fitting procedure 


(IX) 


0.60 
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1.06 
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0.47 


0.42 


0.40 
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0.54 
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(X) 


0.70 
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0.78 
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(XI) 


1.07 


0.27 


0.23 
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0.16 


0.07 


0.09 


0.16 


0.28 


0.16 


0.14 


0.14 


0.53 


(Xll) 


-0.49 


0.01 
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-0.04 


-0.03 


0.12 





0.33 


0.50 


0.10 
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0.14 


0.27 


^syst 


3.70 
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0.94 
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plitude. 



C. Ss ~ Sp measurement 



As explained in previous sections, measurements are 
sensitive to the phase difference between S- and P-waves. 
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FIG. 16: (color online) Projections of background subtracted data (black dots) and fitted MC signal distributions (hashed 
histogram) versus the four kinematic variables in the mass region between 800 and 900 MeV/c'^. Error bars correspond to 
statistical uncertainties only. The signal contains S-wetve, K (892)" and K (1410)° components. Lower plots are the ratio 
between data and the fitted signal. 



TABLE XIII: Positions of the center and values of half the 
mass intervals used in the phase measurement. 

mass bin ttikit (GeV/c^) ArriKTr ( GeV/c^) 



1 


0.707 


0.019 


2 


0.761 


0.035 


3 


0.828 


0.032 


4 


0.880 


0.020 


5 


0.955 


0.055 


6 


1.047 


0.037 


7 


1.125 


0.041 


8 


1.205 


0.039 


9 


1.422 


0.178 



This quantity is given in Fig. 21 for different values of the 
Ktt mass using results from the fit explained in Section 
VIII B. Similar values are obtained if the K (1410)° is 
not included in the P-wave. 



D. Search for a D-wave component 

A D-wave component, assumed to correspond to the 
(1430)", is added in the signal model using expres- 
sions given in Eq. (21) and (30-32). As the phase of the 



TABLE XIV: Fit results for a signal made of S'-wave, 
IT* (892)° and !^*(1410)° components. The S-wave phase is 
measured for several values of the K-k mass and its amplitude 
is parameterized according to Eq. (28) . The two last columns 
give the values of the P-wave phase, which includes K (892)° 
and K (1410)° components, and the values of the difference 
between the S- and P-wave phases. Quoted uncertainties are 
statistical only, systematic uncertainties are given in Table 
XII. The same uncertainties apply to 5$ and 5s — 5p- 



variable 


result 






^A'* (1410)0 


0.079 ± 0.004 






<5if (1410)0 (°) 


-8.9 ±21.5 






rs 


0.463 ± 0.068 








0.21 ±0.18 












Ss-SpC) 


5i 


16.8 ± 11.7 


2.0 


14.8 


52 


31.3 ±5.5 
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47.7 ± 1.4 
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55.0 ±4.2 


161.4 
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71.2 ±6.9 
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FIG. 17: (color online) Projections of background subtracted data (black dots) and fitted MC signal distributions (hashed 
histogram) versus the four kinematic variables in the mass region between 900 and 1000 MeV/c^. Error bars correspond to 
statistical uncertainties only. The signal contains S-wave, K (892)" and K (1410)° components. Lower plots are the ratio 
between data and the fitted signal. 



If* (1410)°, relative to the K*{892)°, is compatible with 
zero this value is imposed in the fit. For the Z?-wave, 
its phase (Sd) is allovifed to be zero or tt. Fit results 
are given in the last column of Table VIII. The total 
value is 2888 and the number of degrees of freedom 
is 2786. This corresponds to a probability of 8.8%. The 
value zero is favored for Sd- The fraction of the decay 
rate which corresponds to the D wave is given in Table 
IX and is similar to the K (1410)" fraction. 



IX. DECAY RATE MEASUREMENT 



The — > A' ^7r+e^i^e(7) branching fraction is mea- 
sured relative to the reference decay channel, 
K~TT^7T~^ {j). Specifically, in Eq. (41) we compare the 
ratio of rates for the decays K^n^e'^i'ei'y) and 

iir~7r+7r+(7) in data and simulated events: this 
way, many systematic uncertainties cancel: 



Rl 



B{D+ ^ K-7V+e+U,)iata 



N{D+ 



> K-n + -K+)data 



(40) 



t(K-1T+e+Ve)data 



C{KTVTT)dat 



N(D+ K-K+TT+)data C{KTiev)data ' 



Introducing the reconstruction efficiency measured for 
the two channels with simulated events, this expression 
can be written: 



Rd 



e I'eldata 



K 7r+7r^ 



C{K'IVK)data 
C{K-KeVf,)data 



e{K ■K+e+Ve)data 



e{K 7r+7r+)AfC 



(41) 



The first line in this expression is the product between 
the ratio of measured number of signal events in data for 
the semileptonic and hadronic channels, and the ratio 
of the corresponding integrated luminosities analyzed for 
the two channels: 



£{KTVeUe)data 



98.7 fb- 
100.5 fb" 



(42) 



The second line of Eq. (41) corresponds to the ratio be- 
tween efficiencies in data and in simulation, for the two 
channels. The last line is the ratio between efficiencies for 
the two channels measured using simulated events. Con- 
sidering that a special event sample is generated for the 
semileptonic decay channel, in which each event contains 



a decay Z?^ 



+ 



K 



■*o 



K 7r+, whereas the 
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FIG. 18: (color online) Projections of background subtracted data (black dots) and fitted MC signal distributions (hashed 
histogram) versus the four kinematic variables in the mass region between 1000 and 1600 MeV/c'^. Error bars correspond to 
statistical uncertainties only. The signal contains S'-wave, K (892)" and K (1410)° components. Lower plots are the ratio 
between data and the fitted signal. 



D'^ — )■ K TT^TT"*" is reconstructed using the e^e cc 
generic simulation, the last term in Eq. (41) is written: 

€{K-7t+7v+)mc ^ N{D+ ^ K-n+n+)MC ^^3^ 
e(K*°e+iy,)MC N{D+ ^W°e+iy,)MC 

^ jV(D+->Z*°e+^.)i- 

2N{cc)k^^V{c ^ D+)B{D+ ^ K-TT+n+)MC 

where: 

• N{D+ K*°e+Ve)tic = 1-17 X is the number 
of generated signal events; 

• N{cc)k-ktt = 1.517 X 10^ is the number of 
e^e~ — ?> cc events analyzed to reconstruct the 
D+ —J' K~Tr^'K'^ channel; 

• V{c D+) = 26.0% is the probability that a 
c-quark hadronizes into a in simulated events. 
The is prompt or is cascading from a higher 
mass charm resonance; 

• B{D+ K-Tr+TT+)Mc = 0.0923 is the branching 
fraction used in the simulation. 



A. Selection of candidate signal events 

To minimize systematic uncertainties, common selec- 
tion criteria are used, as much as possible, to reconstruct 
the two decay channels. 



1. The D — > K TT TV decay channel 

As compared to the semilcptonic decay channel, the se- 
lection criteria described in Section VI A are used, apart 
for those involving the lepton. The number of signal 
candidates is measured from the K~T:^Tr~^ mass dis- 
tribution, after subtraction of events situated in side- 
bands. The signal region corresponds to the mass in- 
terval [1.849, 1.889] GeV/c^ whereas sidebands are se- 
lected within [1.798, 1.838] and [1.900, 1.940] GeV/c^. 
Results are given in Table XV and an example of the 
K~tt^t:^ mass distribution measured on data is dis- 
played in Fig. 22. 

The following differences between data and simulation 
are considered: 

• the signal mass interval. 

Procedures have been defined in Section VI C 5 such 
that the average mass and width of the D+ — )■ 
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FIG. 19: (color online) Points (full circles) give the S-wave phase variation assuming a signal containing S-wave, K (892)° 
and K (1410)° components. The S-wave phase is assumed to be constant within each considered mass interval and parameters 
of the K (892)° are fixed to the values given in the third column of Tab. VIII. Error bars include systematic uncertainties. 
The full line corresponds to the parameterized 5-wave phase variation obtained from the values of the parameters quoted in 
the same column of Table VIII. The phase variation measured in Kit scattering by Ref. [4] (triangles) and LASS [5] (squares), 
after correcting for 5^''^, are given. 



TABLE XV: Measured numbers of signal events in data and 
simulation satisfying F^c > 0.5. 

Channel Data Simulation 
K'Ti+e+u^ 70549 ±363 330969 
K~-K+-K+ 52595 ± 251 68468 ± 283 



tracking efficiency cancel when considering the si- 
multaneous reconstruction of the pion and the elec- 
tron. What remains is the difference between data 
and simulation for electron identification which is 
included in the evaluation of systematic uncertain- 
ties. 



K TT'^TT'^ reconstructed signal in data and simu- 
lated events are similar. 

• the Dalitz plot model. Simulated events are gener- 
ated using a model which differs from present mea- 
surements of the event distribution over the Dalitz 
plane. Measurements from CLEO-c [9] are used 
to reweight simulated events and we measure that 
the number of reconstructed signal events changes 
by a factor 1.0017 ± 0.0038. This small variation 
is due to the approximately uniform acceptance of 
the analysis for this channel. 

• the pion track. As compared with the K~''r:'^ e'^Ve 
final state, there is a 7r+ in place of the e+ in the 
reference channel. As there is no requirement on 
the PID for this pion we have considered that pos- 
sible differences between data and simulation on 



2. The -D — > K -R^e decay channel 

The same data sample as used to measure the D'^ — >■ 
i^T" 7r+7r+ is analyzed. Signal events are fitted as in Sec- 
tion VII. The stability of the measurement is verified 
versus the value of the cut on Fee which is varied be- 
tween 0.4 and 0.7. Over this range the number of signal 
and background events change by factors 0.62 and 0.36 
respectively. The variation of the ratio between the num- 
ber of selected events, 

^ _ N{D+ ^ K-TT+^+)mC N{D+ ^ K-TT+e+V,)data 
N{D+ ^K*\+Ve)MC N{D+ K-TT+TT+)data 

(44) 

in data and simulation is given in Table XVI. 

Relative to the value for the nominal cut {Fee > 0-5), 
the value of Rn for Fee > 0.4 is higher by 0.00038 ± 
0.00063 and for Fee > 0.7 it is higher by 0.0018 ± 0.0011. 
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FIG. 20: (color online) Comparison between present measurements of the 1=1/2 S— wave phase variation with the Ktt mass 
and previous resuhs from Estabrooks et ah [4], LASS [5], E791 [6], FOCUS [7, 8], and CLEO [9]. 



TABLE XVI: Variation of the ratio between the numbers of selected events in data and simulation for different values of the 
cut on Fee- 









Fee > 0.4 


Fee > 0.5 


Fee > 0.7 


N{D+ - 


K--JT + TV + 


'MC 


72206 ± 292 


60468 ± 283 


59259 ± 259 


N{D+ - 


-> K-TY + TV+] 


data 


55361 ± 260 


52595 ± 251 


45627 ± 230 


N{D+ 




MC 


381707 


330969 


237104 


N{D+ - 




) data 


81322 ± 383 


70549 ± 363 


50989 ± 303 




Rn 




0.2779 ± 0.022 


0.2775 ± 0.0023 


0.2793 ± 0.0026 



Quoted uncertainties take into account events that are 
common when comparing the samples. These variations 
are compatible with statistical fluctuations and no addi- 
tional systematic uncertainty is included. 

To select semileptonic decay candidates a cut is applied 
on the probability of the mass-constrained fit at 0.01. 



In a previous analysis of the decay — > K^e^h'e [41] we 
measured a value of 1.0062±0.0006 for the ratio between 
the efficiency of this cut in simulation and data. We use 
the same value in the present analysis because this prob- 
ability depends on the capability to reconstruct the D 
direction and momentum and to estimate corresponding 
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FIG. 21: (color online) Difference between the 7 = 1/2 5'- and P-wave phase versus the Kn mass. Measurements are similar 
whether or not the K (1410)'' is included in the P-wave parameterization. Results are compared with measurements from K-k 
scattering [4, 5]. The continuous and dashed lines give the phase variation with a minus sign {—5p) for the K (892)" and 
K (892)" + K (1410)", respectively. The difference between these curves and the measured points corresponds to the S'-wave 
contribution. 



event. 




1.85 1.9 

mta(GeV/c') 

FIG. 22: K~-K^Tx^ mass distribution measured in data. The 
signal and sideband regions are indicated. 



uncertainties on these quantities which are obtained, not 
from the studied decay channel, but from the rest of the 



B. Decay rate measurement 

Measurement of the _D+ X^7r+e+j/e branching 
fraction and of the contributing 5— wave, K (892)", 
K (1410)" and i^2(1430)" components is important to 
verify if the sum of exclusive channels in Z3-meson 
semileptonic decays agrees with the inclusive value. From 
the measurement oiB{D^ — > K (892)"e"'"i'e) the value of 
1^1 (0)1 is obtained and provides, with r2 and ry, the ab- 
solute normalization for the corresponding hadronic form 
factors. These values can be compared with Lattice QCD 
determinations. 

Combining all measured quantities in Eq. (41), the rel- 
ative decay rate is: 



Rd = 0.4380 ± 0.0036 ± 0.0042 



(45) 



where uncertainties are statistical and systematic respec- 
tively. Using the CLEO-c value for the branching fraction 
B{D+ K-Tr+TT+) = (9.14 ± 0.20)% [52], gives: 

B{D+ R-TT+e+Ve) = (4.00±0.03±0.04±0.09) x 10"^ 

(46) 

where the last quoted uncertainty comes from the ac- 
curacy of B{D~^ — )■ if~7r+7r+). To evaluate the con- 
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F*0 



tribution from the K " , results obtained with the S + 
IT* (892)0 _^ 7?* (1410)° signal model are used. The 
branching fraction for £)+ — ^ K e+z^e is obtained after 
subtracting the S- and K (1410)°-wave contributions: 



B{D+ K*"e+iye) x B{K*" K-tt+) 
(3.77 ± 0.04 ± 0.05 ± 0.09) x lO'^. 



(47) 



The last uncertainty corresponds to external inputs. 

The corresponding value of ^i(O) is obtained by inte- 
grating Eq. (16), restricted to the K contribution, over 
the three angles: 



dr 



1G| 



dq^dm^ 3 (47r 



,5 2 



D 



3{i^iip+i^2ii^+i^3ir} 



(48) 



Assuming that the if* (892)° meson is infinitely nar- 
row, integrating over the remaining variables gives: 



r ^ nB{D+ ^K*%+i.,)B(K*° ^ K-n+) ^^^^ 



G|||K.II'2 



td+ 

2-1 



with 



-I 



Q^max 2 



r+\ 



X \\Ho\^ + \H+\^ + \H_\^]dq^ 



(50) 



and: 



Ai (0) = 0.6200 ± 0.0056 ± 0.0065 ± 0.0071. (51) 

For this last evaluation, the values T£)+ = (10.40 ± 
0.07) X 10-"s for the D+ lifetime [16] and \Vcs\ = 
0.9729 ± 0.0003 are used. Corresponding uncertainties 
are included in the last quoted error in Eq. (51). 

If instead of considering a if* (892)° with zero width, 
the fitted mass distribution of the resonance is used in 
the integral of the differential decay rate versus q'^ and 
m^, the form factor normalization becomes: 

^i(0)|g2,„2 = 0.9174 ±0.0084 ±0.0097 ±0.0105. (52) 

This value depends also on the normalization adopted for 
the mass distribution which is given in Eq. (26). 



— *o 

system in this decay is dominated by the K meson. In 
— *o 

addition to the K meson we measure a contribution of 
the K-1T+ S-wave component of (5.79 ± 0.16 ± 0.15)%. 
We find a small contribution from the K (1410)° equal 
to (0.33±0.13±0.19)%. This value agrees with the naive 
expectation based on corresponding measurements in r 
decays. The relative phase between the K (892)° and 
K (1410)° components is compatible with zero whereas 
there is a negative sign between the 5- and P-wave am- 
plitudes. A fit to data of similar probability is obtained 
including a £>-wave component with a fraction equal to 
(0.19±0.09±0.09)%. In this case the 1^* (1410)° fraction 
becomes (0.16 ±0.08 ±0.14)%. As these two components 
do not exceed a significance of three standard deviations, 
upper limits at the 90% C.L. are quoted in Table XVII. 

Using a model for signal which includes S'-wave, 
K*(892)°- and !K*(1410)° contributions, hadronic form 

factor parameters of the K* component are obtained 
from a fit to the five-dimensional decay distribution, as- 
suming single pole dominance: ry = V{0)/Ai{0) = 
1.463±0.017±0.032, r2 = A2(0)/Ai(0) = 0.801±0.020± 
0.020 and the pole mass of the axial vector form fac- 
tors ruA = (2.63 ±0.10 ±0.13) GeV/c^. For compari- 
son with previous measurements we also perform a fit 
to data with fixed pole mass ttia = 2.5 GeV/c^ and 
mv = 2.1 GeV/c^ and including only the S and K* (892)° 
signal components; it gives ry = 1.493 ± 0.014 ± 0.021 
and r2 = 0.775 ± 0.011 ± 0.011. 

We have measured the phase of the S'-wave compo- 
nent for several values of the ii'~7r"'" mass. Contrary 
to similar analyses using charm meson decays, as in 
£)+ — )• if~7r+7r+, we find agreement with correspond- 
ing S'-wave phase measurements done in K~p interac- 
tions producing K~Tr~^ at small transfer. This is a con- 
firmation of these results and illustrates the importance 
of final state interactions in L'-meson hadronic decays. 
As compared with elastic if "tt"*" scattering there is an 
additional negative sign between the 5- and P-wave, in 
the semileptonic decay channel. This observation 
does not contradict the Watson theorem. We have deter- 
mined the parameters of the if * (892)° meson and found, 
in particular, a width smaller than the value quoted in 
[16]. Our result agrees with recent measurements from 
FOCUS [14], CLEO [9] and, r decays (for the charged 
mode) [12]. Comparison between these measurements 
and present world average values is illustrated in Tab. 
XVII. Our measurements of the S-wave phase have 
large uncertainties in the threshold region and it remains 
to evaluate how they can improve the determination of 
chiral parameters using, for instance, the framework ex- 
plained in ref. [19]. 



X. SUMMARY 

We have studied the decay £>+ — )■ if ~7r+e+i/e with a 
sample of approximately 244 x 10^ signal events, which 
greatly exceeds any previous measurement. The hadronic 
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TABLE XVII: Comparison between these measurements and present world average results. Values for B{D'^ 
K (1410)°/it'2(1430)''e"'"^e) are corrected for their respective branching fractions into K~'k~^ . 



Measured quantity 


This analysis 


World average [16] 


rnjf* (892)0 (MoV/c^) 
Tk. (892)0 (McV/c^) 
rBw{Ge\/c)-^ 


895.4 ± 0.2 ± 0.2 
46.5 ± 0.3 ± 0.2 

2.1 ±0.5 ±0.5 


895.94 ± 0.22 
48.7 ± 0.8 

3.4 ± 0.7 [5] 


rv 
r2 

mA(GeV/c^) 


1.463 ±0.017 ±0.031 
0.801 ± 0.020 ±0.020 
2.63 ±0.10 ±0.13 


1.62 ±0.08 
0.83 ± 0.05 
no result 



B{D+ K-'K+e+v^){%) 4.00 ± 0.03 ± 0.04 ± 0.09 4.1 ± 0.6 

B{D+ if-7r+e+i^e)K* (892)0 (%) 3.77 ± 0.04 ± 0.05 ± 0.09 3.68 ± 0.21 

B{D+ K-'K+e+v^)s_^av^{V(,) 0.232 ± 0.007 ± 0.007 ± 0.005 0.21 ± 0.06 

B{D+ ^*(1410)°e+i^e)(%) 0.30 ± 0.12 ± 0.18 ± 0.06 (< 0.6 at 90% C.L.) 

B{D+ -» ;g2(1430)°e+i/e)(%) 0.023 ± 0.011 ± 0.011 ± 0.001 (< 0.05 at 90% C.L.) 
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Appendix A: Error matrices for the nominal fit 



The correlation matrix between statistical uncertain- 
ties is given in Tab. XVIII for parameters fitted using 
the nominal model. Statistical uncertainties are quoted 
on the diagonal. The elements of the statistical error 
matrix are equal to pij (j,; aj , where pij is an off diagonal 
element or is equal to 1 on the diagonal. 

In Tab. XI, systematic uncertainties quoted in lines 
labelled III, VII and XI are the result of several contri- 
butions, combined in quadrature. In Tab. XIX these 
components are detailed because each contribution can 
induce a positive or a negative variation of the fitted 
quantities. 

The error matrix for systematic uncertainties on fit- 
ted parameters is obtained using values of the variations 
given in Tab. XI and XIX. For each individual source 
of systematic uncertainty we create a matrix of elements 
equal to the prodiict 5i Sj of the variations observed on 
the values of the fitted parameters i and j. For systematic 
uncertainties IX and X, which have a statistical origin, 
we multiply thesc^ quantities by the corresponding ele- 
ments of the statistical correlation matrix (Tab. XVIII). 
These matrices are summed to obtain the matrix given 
in Tab. XX. 
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TABLE XVIII: Correlation matrix for the S + K* (892)° + 7f*(1410)" nominal fit. On the diagonal, statistical uncertainties 
(crl*"* ) of fitted quantities (i) are given. 



AM^, (892)0 Ar^*(g92)0 AfBlV 


AniA 


Arv 


Ar2 


Ars 


Ar« 


A 1/2 

^°-S,BG 


^''if* (1410)0 


ASjf* (1410)" 


ANs 


ANb 


0.211 0.656 -0.842 


-0.158 


0.142 


-0.131 


0.116 


-0.043 


-0.673 


0.899 


-0.774 


-0.254 


0.304 


0.315 -0.614 


-0.002 


0.020 


0.007 


0.007 


0.024 


-0.470 


0.624 


-0.632 


0.027 


-0.021 


0.476 


0.163 


-0.165 


0.141 


-0.347 


0.270 


0.657 


-0.907 


0.846 


0.334 


-0.394 




0.0972 


-0.548 


0.840 


-0.045 


-0.070 


0.065 


-0.080 


0.087 


0.099 


-0.118 






0.0166 


-0.518 


0.048 


0.034 


-0.060 


0.101 


-0.126 


-0.124 


0.136 








0.0201 


-0.016 


-0.080 


0.051 


-0.058 


0.080 


0.116 


-0.133 










0.0286 


-0.968 


-0.157 


0.191 


-0.048 


-0.136 


0.159 












0.0640 


0.130 


-0.133 


-0.043 


0.146 


-0.173 














0.138 


-0.767 


0.396 


0.148 


-0.179 
















0.0163 


-0.721 


-0.269 


0.318 


















13.0 


0.288 


-0.336 




















713.0 


-0.609 






















613.2 



TABLE XIX: Systematic uncertainties on parameters fitted using a model for the signal which contains S'-wave, K (892) 
and K (1410)° components, expressed as (a;[0] — x[i])/aatat- (Ilia) uncertainty associated with electron identification, (Illb) 
uncertainty associated with kaon identification, (Vila) : pole mass changed by -30 MeV/c^ for the decay channel D° — ^ K~e^Ve, 
(Vllb) : pole mass changed by -100 McV/c^ for scmilcptonic decays of and into a pseudoscalar meson, (VIIc) : pole 
mass changed by 100 MeV/c^ for Ds meson scmilcptonic decays, (Vlld-j) refer to D V decays; (Vlld) : r2 changed 

from 0.80 to 0.85, (Vile) : rv changed from 1.50 to 1.55, (Vllf) : mA changed from 2.5 to 2.2 GeV/c^ (Vllg) : rav changed 
from 2.1 to 1.9 GeV/c^ (Vllh) sbw changed from 3.0 to 3.3 GeV-\ (Vlli) : r^.(892)o varied by -0.5 MeV/c^ (XIa-d) : 
it'* (1410)° mass and width, 7^0(1430)° mass and width in this order and using the variations given in Tab. I and, (Xle) : my 
changed by 100 MeV/c^. 



variation AMj^.jggjjo 




Arnw 


AniA 


Arv 


Ar2 


Ars 


Arf 




'"if* (1410)0 


^^Jf* (1410)0 


ANs 


ANb 


III 


a 0.02 
b 0.00 


0.03 
-0.10 


-0.03 
0.05 


0.70 
-0.07 


-0.73 
-0.05 


0.50 
-0.17 


0.12 
-0.07 


-0.07 
0.01 


0.30 
-0.28 


0.07 
0.05 


-0.20 
0.10 


0.07 
-0.07 


-0.08 
0.08 


VII 



a 0.06 0.06 

b 0.01 0.01 

c -0.01 0.00 

d -0.03 0.00 

e -0.03 -0.01 

f 0.00 0.03 

g -0.04 -0.01 

h -0.06 -0.04 

i -0.02 -0.02 



XI 


a 


-0.08 


-0.02 


0.08 


0.00 


0.01 


0.00 


-0.35 


0.34 


-1.00 


-0.34 


-0.20 


-0.04 


0.05 


b 


0.20 


0.09 


-0.21 


0.03 


-0.03 


0.03 


0.44 


-0.41 


3.09 


0.62 


0.28 


0.07 


-0.08 


c 


-0.06 


-0.02 


0.08 


0.01 


-0.01 


0.01 


0.08 


-0.10 


0.05 


-0.39 


-0.08 


0.02 


-0.02 


d 


-0.16 


-0.07 


0.17 


0.01 


-0.01 


0.01 


0.05 


-0.06 


0.11 


0.37 


0.18 


0.02 


-0.03 


e 


0.00 


0.00 


-0.01 


-0.06 


-1.15 


-0.07 


0.04 


-0.03 


-0.01 


0.00 


0.00 


0.01 


-0.01 



0.05 


-0.10 


-0.05 


0.01 


0.17 


-0.15 


0.01 


0.07 


0.02 


0.05 


-0.06 


0.01 


0.00 


-0.01 


0.01 


0.01 


-0.01 


0.00 


0.01 


-0.01 


0.00 


-0.01 


0.01 


-0.01 


0.00 


0.00 


-0.01 


0.01 


0.00 


-0.01 


0.00 


0.00 


0.01 


0.06 


0.00 


0.02 


0.05 


0.01 


-0.02 


0.02 


-0.02 


0.04 


0.06 


-0.07 


0.05 


0.03 


0.07 


0.03 


-0.02 


0.01 


0.01 


-0.03 


0.04 


0.01 


-0.01 


0.07 


-0.16 


0.07 


-0.03 


0.01 


-0.01 


0.05 


0.01 


0.04 


0.12 


-0.15 


0.06 


0.03 


0.08 


0.04 


-0.02 


0.01 


0.02 


-0.04 


0.05 


0.02 


-0.03 


0.02 


0.00 


0.00 


0.00 


0.00 


0.00 


0.02 


-0.04 


0.06 


0.00 


0.00 


0.01 


-0.00 


0.00 


0.01 


0.00 


-0.01 


0.00 


-0.00 


0.00 


0.01 


-0 01 
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TABLE XX: Correlation matrix for systematic uncertainties of the S + K (892)° + if* (1410)° nominal fit. On the diagonal, 
total systematic uncertainties {cr"^"*') on fitted quantities (i), are given. 





AniA 


Arv 




Ars 


Ar« 


A 

S BG 


ATj^* (1410)0 


^^K* (1410)'^ 


ANs 


AA^B 


0.226 0.569 -0.827 
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0.731 
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